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“methodology  is  to  use  the  logic  of  the  computet  to  take  a  scalar  continuous  out¬ 
put  signal,  sample  it  at  a  relatively  high  rate,  and  then  sort  out  the  various 
samples  into  a  vector  of  ^pseudo  measurements  X  The  potential  of  such  a  scheme 
is  significant,  since  it  offers,  among  other  things,  the  possibility  of  using 
fewer  sensors  and  the  notential  for  providing  gracefully  degrading  system  per¬ 
formance  in  the  face  of  sensor  failure. 

_ _ _ J 

The  pseudo  measurement  concept  may  be  likened  to  the  observer  concept  of 
linear  system  theory.  An  "observer"  is  a  dynamic  subsystem  which  processes  the 
available  measurements  and  control  Inputs  to  obtain  asymptotic  estimates  of 
those  states  which  are  not  measured.  The  distinction  between  the  observer  and 
pseudo  measurement  concept  is  as  follows.  Additional  pseudo  states  required  for 
proper  control  action  are  obtained  via  design  of  a  multi-rate  or  multiple-order 
sampling  sequence  in  the  pseudo  mease. ement  approach,  while  the  estimate  of 
system  states  not  measured  is  attained  by  design  of  a  dynamic  subsystem  in  the 
observer  approach.  The  pseudo  measurement  concept  is  based  on  using  information 
contained  in  a  scalar  measurement,  sampled  at  a  high  rate,  to  generate  a  lower 
rate  measurement  vector .  This  vector  is  likened  to  pseudo  states  since  it 
effectively  provides  distinctively  clfferent  information  in  the  form  of  separate 
sampling  sequences  shifted  in  time. 


-^The  report  is  organized  in  three  volumes.  Volume  I  contains  the  theoretical 
developments  as  well  as  illustrative  examples  and  case  studies.  The  basic  tools 
needed  to  assess  the  pseudo  measurement  concept  are.  reviewed  and  extended.  The 
concept  of  the  frequency  response  of  the  continuous  variables  of  single-rate 
digitally  controlled  system  is  timi  given  sufficient  generality  to  treat  multi¬ 
rate  systems  (including  the  pseudo  measurement  concept).  A  YF-16  case  study  is 
used  to  compare  pseudo  measurement  gains  against  conventional  YF-16  analog  con¬ 
troller  gains.  The  pseudo  wear urc  lent  gains  are  found  to  be  comparable,  with  the 
conventional  gains,  indicating  tin-  approach  does  not  require  control  authority 
more  excessive  than  an  analog  controller  design. a  _ 

A* 

As  a  "byproduct,"  an  .'.mi  or  taut  relationship  is'deduced  which  made  it  feasible 
!  to  apply  the  switch  decompe:- ifi.cn/ frequency  response  technique  to  the  error 
analysis  of  simulations  of  continuous  systems  and/or  discretely  controlled  con¬ 
tinuous  systems.  An  illustrative  example  demonstrates  the  significant  spectral 
differences  between  a  closed-loop  system  (employing  an  analog  or  digital  con- 
Liuller)  and  an  all-digital  simulation  of  the  closed-loop  system.  The  example 
also  treats  problems  encountered  when  simulation  software  is  implemented  on  two 
different  computers,  each  working  in  a  slightly  different  frame  time  (the  so- 
called  Independent  processor  problem).  This  provides  a  tool  which  can  be  profit¬ 
ably  applied  to  the  "error"  analysis  of  simulations;  perhaps  to  predict  expected 
differences  between  actual  flight  test  and  man-in-the-loop  moving-base  simula¬ 
tions. 

Volume  II  describes  two  algorithms  useful  in  the  analysis  of  multi-rate 
systems,  while  Volume  III  gives  FORTRAN  listings  for  these  algorithms.  The 
first,  algorithm  converts  a  continuous  transfer  function  in  the  s-domain  into 
alternate  descriptions  in  the  z-,  w-,  and  w'--domains.  The  second  algorithm  j 
calculates  a  low  data  rate  discrete  transform  given  a  high  data  rate  discrete 
transfer  function.  This  algorithm  is  presented  in  Section  IV  of  Volume  I  and  | 
later  derived  in  Section  IV  of  Volume  II, 
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SECTION  I 


INTRO DUCT ION 


Many  current  digital  design  procedures  attempt  only  to  approach  the 
characteristics  of  existing  analog  designs.  No  explicit  methodology 
exists  which  can  exploit  the  unique  characteristics  of  multi-rate  or 
multiple-order  digital  computation  and  filtering.  The  objective  of  this 
research  effort  is  to  develop  a  multi- rate/multiple-order  theory  for  the 
design  of  closed-loop  digital  control  systems  with  characteristics  and 
attributes  that  are  either  difficult  or  impossible  to  attain  with  analog 
or  single-rate  controllers.  The  key  idea  in  this  methodology  is  to  use 
the  logic  of  the  computer  to  take  a  scalar  continuous  output  signal, 
sample  it  at  a  relatively  high  rate,  and  then  sort  out  the  various  sam¬ 
ples  into  a  vector  of  "pseudo  measurements."  The  potential  of  such  a 
scheme  is  significant,  since  it  offers,  among  other  things,  the  possi¬ 
bility  of  using  fewer  sensors  and  the  potential  for  providing  gracefully 
degrading  system  performance  in  the  face  of  sensor  failure. 

The  pseudo  measurement  concept  may  be  .likened  to  the  observer  con¬ 
cept  of  linear  system  theory.  Recall  that  an  "observer"  is  a  dynamic 
subsystem  which  processes  the  available  measurements  and  control  inputs 
to  obtain  asymptotic  estimates  of  those  states  which  are  not  measured. 
The  distinction  between  the  observer  and  pseudo  measurement  concept  is 
as  follows.  Additional  pseudo  states  required  for  proper  control  action 
are  obtained  via  design  of  a  multi-rate  or  multiple-order  sampling  se¬ 
quence  in  the  pseudo  measurement  approach,  while  the  estimate  of  system 
states  not  measured  is  attained  by  design  of:  a  dynamic  subsystem  in  the 
observer  approach.  The  pseudo  measurement  concept  is  based  on  using 
Information  contained  in  a  scalar  measurement,  sampled  at  a  high  rate, 
to  generate  a  lower  race  measurement  vector.  This  vector  is  likened  to 
pseudo  states  since  it  effectively  provides  distinctively  different 
information  in  the  form  of  separate  sampling  sequences  shifted  in  time. 
However,  it  should  be  noted  that  the  pseudo  measurement  co..  “pt  is  not  a 
state  estimator,  but  a  new  design  concept. 
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Several  circumstances  arise  in  applications  whore  the  pseudo  mea¬ 
surement  concept,  by  exploiting  the  unique  characteristics  of  multi-rate 
and  multiple-order  sampling,  can  have  considerable  favorable  impact. 

•  It  is  desirable  to  minimize  the  number  of  sen¬ 
sors,  or  to 

•  Provide  for  gracefully  degrading  system  perfor¬ 
mance  in  the  face  of  sensor  failure. 

•  Not  all  the  states  are  measurable. 

•  Some  measurements  are  noise  free. 

•  Outer  navigation  loops  are  noise  free. 

•  Outer  navigation  loops  are  updated  at  inherently 
slow  rates. 

When  one  or  more  of  these  circumstances  prevail  in  a  particular 
application,  the  contribution  of  this  theory  will  be  in  terms  of  more 
simple  control  laws,  increased  system  performance,  or  more  system  flexi¬ 
bility  in  accommodating  failures. 

The  report  is  organized  in  three  volumes.  Volume  1  contains  the 
theoretical  developments  as  well  as  illustrative  examples  and  case 
studies . 

Volume  II  describes  two  algorithms  useful  in  the  analysis  of  multi¬ 
rate  systems,  while  Volume  III  gives  FORTRAN  listings  for  these  algo¬ 
rithms.  The  first  algorithm  converts  a  continuous  transfer  function  in 
the  s-domain  into  alternate  desc  Options  in  the  z- ,  w-,  and  w'-domains. 
The  second  algorithm  calculates  a  low  data  rate  discrete  transform  given 
a  high  data  rate  discrete  transfer  function.  This  algorithm  is 
presented  in  Section  I V  of  Volume  I  and  later  derived  in  Section  IV  of 
Volume  II. 

Volume  I  is  organized  as  follows.  Vector  switch  decomposition  tech¬ 
niques,  which  provide  a  straightforward  model  of  the  pseudo  measurement 
vector,  are  reviewed  and  extended  in  Section  II.  Also,  in  this  section 
the  pseudo  measurement  concept  is  defined  and  illustrated  with  an  exam¬ 
ple  . 
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An  important  facet  of  this  effort  hinges  on  the  ability  to  interpret 
the  properties  of  muLci-rate/multiple-order  closed-loop  systems.  Fre¬ 
quency  response  techniques  provide  effective  tools  to  investigate  the 
spectral  characteristics  of  tuese  systems.  Toward  this  end,  the  fre¬ 
quency  response  of  the  continuous  variables  of  eingle-rate  systems  are 
reviewed  in  Section  III,  with  the  necessary  extensions  for  the  analysis 
of  multi-rate  systems,  given  in  Section  IV. 

The  extensions  given  in  Sectioi  IV  make  it  possible  to  investigate 
the  spectral  content  of  the  continuous  variables  of  multi-rate  and/or 
multiple-order  configured  systems.  In  addition,  an  important  relation¬ 
ship  is  deduced  which  makes  it  feasible  to  apply  the  switch  decoraposi- 
t ion/ frequency  response  technique  to  the  error  analysis  of  simulations 
of  continuous  systems  and/or  discretely  controlled  continuous  systems. 
An  illustrative  example  demonstrates  significant  spectral  differences 
between  a  closed-loop  system  (employing  an  analog  or  digital  controller) 
and  an  all-digital  simulation  of  the  closed-loop  system.  The  example 
also  treats  the  problems  encountered  when  simulation  software  is  imple¬ 
mented  on  two  different  (independent)  computers,  each  working  in  a 
slightly  different  frame  time. 

The  application  of  the  pseudo  measurement  concept  is  discussed 
further  in  Section  V.  In  particular,  a  cast  study  using  a  longitudinal 
model  of  the  Y-16  is  used  to  study  the  practicality  of  the  approach 
(e.g.,  investigate  the  magnitude  of  the  required  feedback  gains). 

A  summary  and  conclusions  arc  presented  in  Section  VI. 
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SECTION  II 


VECTOR  SWITCH  DECOMPOSITION:  CONSTRUCTION 
OF  THE  PSEUDO  MEASUREMENT  VECTOR 


A.  INTRODUCTION 


The  prime  objective  of  this  research  effort  is  the  realization  of 
closed-loop  digitally  controlled  systems  with  desirable  attributes  not 
available  in  continuously  controlled  (analog)  designs.  To  appreciate 
the  considerable  merit  of  this  viewpoint  one  need  only  recall  that 
couplers  (data  holds)  which  tie  the  digital  computer  to  the  control 
actuators  are  themselves  filters  possessing  characteristics  which  are 
unrealizable  with  continuous  linear  filter  sections.  For  example,  the 
slewer  data  hold  forces  the  control  actuator  to  move  with  smooth  ramp- 
like  deflections  which  closely  approximate  the  "continuous"  deflec¬ 
tions.  Yet,  it  inserts  deep,  wide  notches  into  the  frequency  response, 
at  multiples  of  the  sampling  frequency,  which  have  no  counterpart  in 
continuous  linear  filter  theory.  The  impact  of  this  frequency  response 
characteristic  upon  system  performance  can  be  favorable  or  unfavor¬ 
able.  It  can  certainly  have  an  unfavorable  effect  if  unanticipated 
inputs  and/or  control  inputs  are  attempted  or  required  at  o\.  .  of  the 
notch  frequencies.  However,  when  anticipated,  this  frequency  response 
characteristic  may  be  used  to  advantage  (for  suppressing  feedback  of 
structural  modes,  for  example)  and  its  disadvantages  avoided. 

In  short,  the  ability  of  the  digital  computer  to  perform  a  variety 
of  unusual  operations  on  input  data  lends  credence  to  the  idea  of 
achieving  closed- loop  systems  with  unusual  attributes.  This  is  the  key 
point  of  the  present  study  where  the  logic  of  the  computer  is  used  to 
take  a  scalar  continuous  output  signal,  sample  it  at  a  relatively  high 
rate,  and  then  sort  it  out  ( in  a  particular  Eashion)  Into  a  vector  of 
"pseudo"  measurements.  The  potential  of  the  idea  is  significant,  since 
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it  offers  an  alternative  patli  approach  for  providing  gracefully  degrad¬ 
ing  system  performance  in  the  face  of  sensor  failure. 

To  begin,  there  is  a  brief  review  of  the  multi-rate  transform  domain 
(primarily  in  order  to  review  notation).  Then,  to  put  the  pseudo  meas¬ 
urement  concept  in  proper  perspective,  we  discuss  vector  switch  decompo¬ 
sition  anu  show  tiow  it  is  used  to  model  multi-rate  systems.  Following 
this,  a  multiple-order  sampling  model  is  discussed. 

Since  any  viable  digital  theory  must  be  capable  of  modeling  computa¬ 
tional  delay,  a  discussion  of  non-synchronous  sampling  is  included  since 
it  provides  a  basic  theoretical  structure  for  treating  time  delays. 

Completing  this  background  description,  we  are  then  in  a  good  posi¬ 
tion  to  set  forth  the  pseudo  measurement  concept  and  indicate,  via  an 
example,  its  potential  for  achieving  closed-loop  configurations  with 
unusual  attributes. 

B.  REVIEW  OF  THE  MULTi-RAXE  TRANSFORM  DOMAIN 

A  basic  property  of  the  transform  domain  algebra  developed  in  Ref.  1 
is  that  a  high  sampling  rate  "operates  through"  a  low  one,  provided  that 
the  ratio  of  the  higher  to  the  lower  is  an  integer  value. 


R  /  RT/M 

C  /  c™ 

G 

T/M 

T/N 

Figure  1.  Multi-Rate  System 
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tn  FI,?,.  I,  the  output  equation  is 


C  =  GRT/M  (l) 

Then 

<;T/N  =  [  GR'T/M  1T/N 

=  qT/M  |^T / M  (2) 

where  N/M  is  an  Integer.  The  superscript  notation  denotes  the  sampling 
rate  involved.  For  example,  indicates  that  the  signal  C  is  sampled 

at  N/T  samples  per  second. 

The  general  notation  Implies  either  C^^(s)  or  in  this 

report,  where  represents  the  Laplace  transform  of  the  sampled 

continuous  time  function  c^^Ct)  and  C^^(z)  its  corresponding  7,  trans¬ 
form.  Tire  correspondence  between  a  sampled  function  in  the  s-domain  and 
its  z- transform  is  given  by 

c'r/N(z)  =  I  CT/N(S)  !z=esT/N  O) 

For  continuous,  stationary  linear  systems  CT/,^(s)  is  in  general  a  non- 
algcbraic  function  of  s  and  Cr^N(z)  is  obtained  from  Cl/M(s)  using  the 
simple  change  of  variable  z  =  e3^^.  However,  could  also  repre¬ 

sent  the  z-domain  model  of  a  completely  discrete  function  (e.g.,  digital 
computations  in  a  computer)  and  thus  have  no  direct  counterpart  in  the 
continuous  domain.  In  analysis  and  design  no  distinction  is  required 
between  the  z-transform  function  derived  from  a  sampled  continuous  func¬ 
tion  [e.g.,  C^^(s)]  and  the  z-’transf orra  function  that  models  a  com¬ 

pletely  discrete  function.  The  z-transform  (and  its  extensions  the  w- 
and  w' -transforms)  thus  provides  a  unified  analysis  and  synthesis  tech¬ 
nique  for  systems  containing  both  continuous  and  dLscrete  elements 
(i.e.,  hybrid  systems).  This  ability  to  model  continuous  and  discrete 
elements  In  a  common  discrete  domain  is  perhaps  the  most  fundamentally 
useful  property  of  the  z-,  w~,  or  w' -domains. 
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The  transform  product  in  Eq .  2  is  obtaired  by  establishing  a  common 
definition  of  the  z-transform  variable.  That  is, 


where 


gT/n(z)rt/m(  )  _  GT/n(Z)rT/M(zN/M) 


=  pS T/M 


The  procedure  is  to  first  calculate  the  T/N  and  T/M  z-transf orms  for 
G(s)  and  R(s),  respectively,  and  then  replace  the  transform  variable 
in  KT/M(zm)  with  zm  =  zN/M. 

Eor  a  given  G  and  R,  a  recursion  equation  can  be  written  for  Eq .  2 
in  a  manner  that  is  similar  to  the  single-rate  case  for  which  N  =  M  •  i. . 
To  see  this,  consider  a  simple  example,  for  which  we  define  z  =  e8^' 


C(s)  = 


cT/N  =  — Ht^rT/m  (?> 

z  “  e  ^/N 

If  the  running  index  of  the  recursion  equation,  n,  is  tied  to  the  z“* 
operator,  the  time  domain  equation  can  be  written  directly  using  the 
fact  that  the  input  is  set  equal  to  R(s)  whenever  n  is  an  integer  value 
of  M/M.  Therefore, 


C  =  p"T/N  (■  + 

un  e  °n_j  m 


Rn  n/(N/M)  =  Integer 


0  n/(N/'M)  *  Integer 


This  agrees  with  the  single-rate  case  when  N  =  M  =  1. 


cn  =  e  T  Cn-1  +  Rn 


These  results  are  appreciably  different  In  the  presence  of  data 
holds.  However,  the  results  remain  quite  tractable  and  are  detailed  in 
Ref.  1  for  the  zero-order  hold  as  well  as  the  slewer  data  hold.  The 
slewer  data  hold  is  a  coupler  which  produces  a  smooth  ramp  output  over 
the  inter-sample  period  and,  in  addition,  provides  a  continuous  output 
even  at  the  sampling  instants.  Transfer  functions  for  representative 
data  holds  are  given  in  Table  1. 


TABLE  1.  SOME  REPRESENTATIVE  DATA  HOLDS 


DATA  HOLD 

- 

TRANSFER  FUNCTION 

Zero-Order 

Hold 

1  -  e~sT 

M0  =  s£ 

First-Order 

Hold 

Ml  =  M^s  + 

Second-Order 

Hold 

M2  =  M(j(s2  +  |f  s  H  p-) 

Triangular 

Data  Hold 

Mo 

Ma  =  y~  esT 

Slewer  Data 

Hold 

m2 

Mg  lew  -p 

The  "T/N"  approach  is  very  useful  in  the  analysis  of  closed-loop 
multi-rate  systems,  since  the  transform  domain  equations  can  be  written 
directly,  without  the  need  for  computing  an  inverse.  Using  these  equa¬ 
tions,  the  time  domain  equations  can  be  written  in  the  same  manner  as  in 
the  single-rate  case.  An  example,  is  described  in  Fig*  2. 


Figure  2.  Closed-Loop  Multi-Rate  System 


First,  writing  the  algebraic  relationships  for  Xj  and  X2  gives 


Xi  = 


s  +  1  1  -  e-sT/3 


s^(s  +  2) 


ET/3 


(10) 


*2 


r 1  _  ~-sT/2i2 


S  i 


(T/2)s2  s2  +  2s  +  5 


X[/ 2 


Ol) 


Recognizing  that 


gT/3  =  rT/3  _  xT/3 


(12) 


and  taking  the  T/N  transform  of  both  sides  produces 


XT/N 


3  ±  1 
s2(s  +  2) 


(RT/3  _  XT/ 3) 


T/N 


(13) 


XT/N 


|  (1  - e~fiT/2) 2 

(  (T/2)s2 


s  +  3 


+  2s  +  5 


T/N 


(14) 
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Since  the  outer  sampling  operator  {*}T/^N  operates  through  the  Inner  one 
If  the  ratio  of  the  inner  sampling  period  to  the  outer  sampling  period, 
T/N,  is  an  integer,  Eqs.  13  and  14  reduce  to 


XT/N 


(  1  -  e  3)  (s  +  1  ) 


s^(s  +  2) 


T/N 


(RT/3  -  XJ/3) 


(15) 


Xj/N 


(l  -  e  sT/2)2(s  +  3) 
(T/2)s2 [s 2  +  2s  +  5) 


XT/2 


(16) 


Observe  that  two  transform  domain  equations  completely  define  the 
closed-loop  fifth-order  system  (exclusive  of  the  data  holds).  Equations 
15  and  16  can  be  used  to  calculate  the  intar-sample  response  with  any 
desired  degree  of  time  resolution  by  increasing  N.  The  only  restriction 
is  that  N  must  be  divisible  by  2  and  3;  otherwise,  it  is  arbitrary. 

The  T/N  approach  is  developed  in  detail  in  Ref.  1.  It  serves  as  the 
basis  for  extending  the  concept  of  frequency  response,  from  merely  being 
the  magnitude  and  phase  of  the  sine  wave  that  fits  the  sample  points  at 
the  sampling  instants,  to  the  case  of  fitting  N  sine  wave.s  to  the  sample 
points  and  N  -  1  inter-sample  points. 

C.  VECTOR  SWITCH  DECOMPOSITION 

AND  MULTI-RATE  SAMPLING 

In  essence,  switch  decomposition  is  a  procedure  wherein  systems 
having  multiple  sampling  operations  (occurring  at  fixed  but  unequal  sam¬ 
pling  intervals  but  with  a  sampling  pattern  which  is  repeated  over  a 
fixed,  finite  time  interval)  are  converted  into  an  equivalent  single 
sample  rate  format.  As  originally  introduced  by  Kranc  (Ref.  2),  the 
method  used  a  summing  point  methodology  which  proved  to  be  extremely 
cumbersome  when  the  ratios  of  the  sampling  periods  become  high.  For 
this  reason  (and  also  because  evolving  state  transition  methods  were 
tending  to  push  transform  methods  into  the  background),  the  method  fell 
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into  disuse.  However,  there  is  much  to  recommend  the  switch  decomposi¬ 
tion  concept  for  use  in  both  time  domain  and  transform  domain  analy¬ 
ses.  In  the  subsection  that  follows  we  will  review  the  basic  concept 
and  remove  some  earlier  restrictions  by  recasting  it  in  /ector  form. 
The  vector  form  simplifies  matrix  block  diagram  manipulation  for  multi- 
loop,  multi-rate  sampled  systems. 

An  example  will  make  the  basic  idea  quite  transparent.  Consider  Lhe 
continuous  signal  shown  in  Fig.  3a  to  be  sampled  at  3/T  samples  per 
second.  This  results  in  the  sample  sequence  shown  in  Fig.  3b.  The 
sampled  values  have  been  numbered  for  easy  reference.  Suppose  we  now 
sample  the  continuous  signal  with  a  sampling  period,  T.  This  results  in 


Figure  3.  Decomposition  of  a  Sample  Sequence 
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•  1  ■  M**»*t*» 


the  sample  sequence  consisting  of  sample  numbers  1,  4,  7,  10,  13,  ... 

shown  in  Fig.  3c.  Define  this  sample  sequence  to  be  R  . 

Next,  advance  the  continuous  signal  R  by  T/3.  Then  sample  the  ad¬ 
vanced  signal  with  a  sampling  period,  T.  This  results  in  a  sample 
sequence  consisting  of  sample  numbers  2,  5,  8,  11,  14,  ...  shown  in  Fig. 
3d.  Define  this  sample  sequence  to  be  (e^^R)^.  Finally,  advance  the 
continuous  signal  R  by  2T/3,  and  sample  it  with  a  sampling  period,  T. 
This  results  In  the  sequence  consisting  of  sample  numbers  3,  6,  9,  12, 
15,  ...  shown  in  Fig.  3e.  Define  this  signal  sequence  to  be  (e2sT/3R)T« 

The  significance  of  the  switch  decomposition  concept  resides  in  its 
ability  to  provide  an  alternative  expression  for  the  original  sequence 

terms  of  several  quantities  which  are  each  sampled  simultane- 

T/  3 

ously  every  T  seconds.  This  alternative  expression  for  R  consists  of 
the  sum  of  RT,  (esT//^R)T,  and 

RT/3  »  RT  +.  (esT/3  R)T  e-sT/3  +  (e2sT/3  R)T  e-2sT/3  (17) 

Equation  17  has  a  simple  factored  equivalent  which  is  the  product  of  two 
vectors  and  the  scalar  R, 


Or,  more  compactly, 

rT/3  ,  W(W*R)t  (19) 

where  . 

W  -  1,  e-sT/3,  e-ZsT/3  (20) 
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and 


1 

esT/3 

e2sT/3 


(21) 


The  extension  of  this  scalar  result  to  the  vector  case  is  discussed 
in  detail  in  Ref.  1  and  will  not  be  covered  here.  We  proceed  to  a  dis¬ 
cussion  of  multiple-order  sampling. 

D.  MULTIPLE-ORDER  SAMPLING 


Suppose  that  a  group  of  sampling  operations  in  a  system  is  repeated 
every  T  seconds,  the  intervals  between  successive  samples  in  the  group 
sequence  being  unequal.  This  is  depicted  in  Fig.  4  for  an  example  case 
of  four  unequally  spaced  samples  occurring  before  the  sampling  sequence 
repeats.  In  Fig.  4  the  notation  X*  is  used  to  denote  that  X  is  sampled 
in  a  multiple-order  manner. 


The  multiple-order  sampling  operation  can  be  modeled  in  terms  of  a 
single-rate  sampling  operation,  using  the  same  type  of  notation  given  in 
Eq.  19.  Let 


X*  =  W(W*X)T 


(22) 
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where 


W  =  [  L,  e“sTl,  e~sT2,  e~sT3 


ii  r  i  sT  ■>  sTn  1 ' 

w*  =  l  1,  e  x,  e  e  J  I 


The  sequence  specified  by  Eq.  22  is  shown  in  Eig.  5.  Note  that 

r  wT  1 


X*  =  W(W  X)T  =  w 


(esTl  X)T 
(esT2  X)T 
(esT3  X)T 


=  XT  +  e_sTl(esTl  X)T  +  e"sT2(esT2  x)T  +  e~sT3(esT3  X)T  (26) 

Tlie  extension  to  the  case  where  X  is  a  vector  is  straightforward.  For 
example,  if  X  is  a  two-dimensional  vector,  let 


Wa 

0 

0 

Wb 

where  the  vectors  W  and  are  each  analogous  to  the  W  vector  defined 
by  Eq .  23  but  each  is  appropriate  to  the  sampling  sequence  for  the  cor¬ 
responding  element  of  the  X  vector. 
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Figure  5.  Components  of  Decomposed  Signal 


E.  NOH- SYNCHRONOUS  SAMPLING 

Non-synchronous  sampling  is  a  basic  tool  useful  for  modeling  dis¬ 
tributed  computation  architecture,  data  skewness  in  the  A/D  and  D/A  con¬ 
version  processes  as  well  as  the  internal  computational  delay  of  the 
digital  computer.  By  definition,  non-synchronous  sampling  occurs  when 
all  the  systems'  sampling  operations  are  repeated  at  the  same  rate  but 
occur  at  different  instants  of  time  (refer  to  Fig.  6). 

In  Fig.  6  both  continuous  signals,  x^  and  are  sampled  at  1/T  Hz, 

but  the  sampler  is  "out  of  sync"  with  the  sampler  by  TQ  seconds. 

The  sampling  operation  for  xj  is  shown  symbolically  in  Fig.  7a  and  for 
x*  in  Fig.  7b.  (*  notation  on  X£  is  used  here  to  indicate  an  "unconven¬ 

tional"  sampling  operation.) 


15 


Figure  6.  A  Set  of  Non-Synchronously  Sampled  Signals 


Figure  7.  Sampling  Notation 


Figure  7b  models  the  non-synchronous  sampler  with,  a  synchronous 

sampler  by  preceding  the  sampler  with  the  operation  W  followed  by  the 

* 

operator  W.  That  is. 


where 


x*  =  w(WAx2)r 


(28) 


e 


-sT0 


sT 

e 


0 


W 


(29) 


Proceeding  according  to  Eq .  29,  one  advances  *2  by  T  seconds,  sam- 
pies  at  the  1/T  rate,  and  then  delays  (W  ^2)  by  10  seconds  to  obtain 
the  time  sequence  (refer  to  Pig.  8).  Note  how  the  non-synchronous  sam¬ 
pling  operations  on  X2  is  modeled  in  terms  of  a  scalar  factor;  thus  the 
dimension  of  the  equivalent  single-rate  sampled  signal,  (W  *2)^1  is  not 
increased . 

The  model  readily  extends  to  the  case  where  x  is  a  vector. 


Figure  8,  Advance,  Sample,  Delay 

V.  THE  PSEUDO  MEASUREMENT  CONCEPT 

It  has  now  been  shown  that  vector  switch  decomposition  can  be  used 
to  model  multi-rate  and  multiple-order  sampling  operations.  Moreover, 
non-synchronous  sampling  (i.c.,  skewed  data)  can  be  treated,  and  there¬ 
fore,  of  course,  computational  time  delays  can  also  be  included.  We  are 
now  in  a  good  position  to  discuss  the  pseudo  measurement  concept,  using 
Fig.  9  to  indicate  the  several  options  available. 
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Figure  9.  Pseudo  Measurement  Block  Diagram 

First  of  all,  the  notation  W^,  W  in  Fig.  9  is  used  to  denote  either 
multi-rite  or  multiple-order  sampling  operations.  The  essence  of  the 
pseudo  measurement  idea  lies  in  the  following  observations  : 

•  X  is  the  multi-rate  (multiple-order)  sampled 
equivalent  of  X.  It  has  the  same  dimension  as 
the  vector  X. 

T 

•  Y  has  tne  row  dimension  of  W  which  may  be  con¬ 
siderably  greater  (but  never  less  than)  the 
dimension  of  the  vector  X. 

Rather  than  feed  back  the  entire  fixed  sequence  of  samples  in  a 
(sealer)  element  of  X* ,  one  may  elect  to  feed  back  selected,  Individual 
samples  In  the  fixed  sequence.  These  individual  samples  are  elements  of 
the  vector  signal,  Y  .  Of  course,  on*  must  exercise  due  regard  for 
physical  realizability  (a  sample  cannot  be  used  before  it  is  availa¬ 
ble).  Nevertheless,  the  digital  computer  has  the  ability  to  sort  out 
individual  samples  from  the  sequence  and  deposit  them  in  their  respec¬ 
tive  appropriate  registers.  This  sorting  process  in  effect  creates  an 
enlarged  vector  of  measurements  at  a  slower  effective  sampling  rate.  We 
shall  call  the  elements  of  this  enlarged  measurement  vector  pseudo 
measurements  because  they  are  merely  alternative  samples  of  the  same 
basic  variables. 

T 

Af  shown  in  Fig.  9,  the  Y  vector  is  defined  in  terms  of  time  ad¬ 
vance  components  and  therefore  caution  must  be  used  in  interpreting  the 
significance  of  the  various  components  (again,  a  sample  cannot  be  used 
before  it  Is  available).  To  avoid  this  difficulty,  the  pseudo  meas¬ 
urement  operation  Is  redefined  in  terms  of  "hat"  operators  (refer  to 
Fig.  10). 
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Figure  10.  Reformulated  Block  Diagram  !' 


If,  for  example, 


=  [  1,  esT/^,  e2sT/3  j' 

and 

W  **  [  1,  e'sT/3j  e"sl’/3  ] 

Then 


(33) 

(34) 


(35) 


Clearly, 


10  0  XT 

X*  =  XT/3  -  W(W*X)T  =  0  e-sT/2  0  (esT/3X)T 

0  0  e-2gT/3  (e2sT/3X)T 


XT 


e-sT/3(esT/3x)T 

e-2sT/3(e2sT/3xjT 


(36) 


Thus,  X  Is  the  desired  vector ;  each  component  has  the  samples  In  the 
correct  temporal  sequence  (refer  to  Fig.  11). 


20 


x> 


X 


,  T/3 


.  T 

x,  =  xT 


£  -«T/3(  iT/3  T 

Xj>  -  e  ( €r  X  ) 


a  =  e'ZsT/3(e2sT/3X  )T 


4 

-X 


V 


1: _ l 


t 


9 

-X 


Figure  11.  X,  X and  the  Components  of  )0 / $ 

C.  AN  ILLUSTRATIVE  EXAMPLE 


Consider  the  single-rate  system  shown  in  Fig.  where  M  is  a  zero- 
order  hold  and  the  plant  is 


3g3  +  15s2  4-  35s  +  13 _ 

(s  +  Dl (s  +  iy-  +  c  3  > 2  ] 


(37) 


Suppose  C(s)  (a  scalar)  is  the  only  output  of  the  system  available  for 
feedback.  Further,  suppose  the  design  objectives  are  to  force  a  unity 
steady-state  response  in  C  to  a  unit  step  input  and  to  obtain  an  as¬ 
signed  set  of  poles  for  the  closed-loop  system.  Let  D(s)  represent 


Figure  12.  A  Single-Rate  Controller,  T  =  1,0  sec 


the  plant  poles,  l.e.,  the  poles  of  G(s).  The  poles  for  the  discretized 
plant  plus  data  hold  are  represented  by  D(z)  where 


<CM)T 


N(z) 

D(z) 


(38) 


and 


D( z)  -  (z  -  0.3679 )[(z  +  0.1340)2  +  (0.019)2]  ,  T  =  1.0  (39) 


with 

[s  +  a]7  H > 


c 


-aT 


[ 


+  a)2  +  b2] 


[  (z  -  e  cos  bT)2  +  (e-a^  sin  bT)2] 
z 2  -  (2e-a-  cos  bT)z  +  e_2a* 


(40) 


';'he  assigned  poles  for  the  closed-loop  system  are  represented  by  A( z) 
where 


A (z)  -  (z  -  0. I353)[ (z  -  0.1988)2  +  (0.3096)2]  (41) 


The  roots  of  A(z)  may  be  thought  of  as  counterparts  to 


A ( s )  =  (s  +  2)[(s  +  l)2  +  (l)2]  (42) 

for  a  continuous  system.  This  particular  assignment  is  quite  arbitrary 
on  our  part  and,  moreover,  we  are  not  particularly  advocating  pole 
assignment  as  a  design  method;  our  objective  is  to  demonstrate  that  the 
design  goals  can  be  achieved  via  the  properties  of  multiple-order  sam¬ 
pling. 

It  may  be  difficult  to  achieve  these  goals  through  design  of  the 
filters  Fi  and  F^.  For  example,  from  the  equation  (see  Fig.  12) 

CT  =  [l  +  (CM)tf{fT]-1  (gm)T  FTRT  (43) 

It  can  be  appreciated  that  the  gain  must  be  split  between  and  F2 

filters  in  order  to  achieve  the  unity  step  response  objective.  More¬ 
over,  it  is  not  readily  apparent  that  one  can  obtain 

D(z)  [  1  +  (CM) r  (F]'F|J  ]  =  A  (z )  (44) 

and  stilJ  maintain  a  third-order  characteristic  equaHon  (an  implicit 
requirement  of  the  pole  assignment  objective). 

Having  considered  a  single-rate  approach,  consider  next  what  multi¬ 
ple-order  sampling  has  to  offer.  A  multiple-order  sampled  configuration 
is  shown  in  Fig.  13.  In  this  simple  configuration,  multiple-order  sam¬ 
pling  is  used  only  on  the  plant  output  variable.  Multiple-order  sam¬ 
pling  could  also  be  employed  for  sampling  the  input,  R,  and  at  the  input 
to  the  data  hold,  U.  Restricting  the  example  to  sampling  the  output 
variable  will  enable  us  to  demonstrate  the  possibilities  of  multiple- 
order  sampling  without  adding  unnecessary  complexity.  The  multiple- 
ord.r  sampling  assumption  allows  us  to  feed  back  the  pseudo  measurement 
vector  b*  instead  of  the  scalar  measurement  C1.  Furthermore,  the  intent 
is  to  use  only  a  ga t n  matrix,  K,  to  attain  the  design  objectives.  The 
particular  choLce  for  W  and  K  is 
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Figure  13.  A  Multiple-Order  Sampling  Controller, 
T  =  1 .0  sec. 


W  =  [  1,  e“(" /10)s>  e~,/is ,  e”  • 6s ,  e'*7s  ]  (45) 

and 

K  =  [  Kj ,  K2,  K3,  K4,  K5  ]  (46) 

Then 
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and 


C 


N 

n 


3s3  +  15s2  4  35b  +  13 
(s  +  l)((s  +  2)2  +  (3)2] 


(59) 


Carrying  out  the  transform  operations  indicated  in  Eq .  57  produces 


Performing  the  matrix  multiplication  and  taking  the  ^averse  in  Eq.  00 
gives 


z  -  i\  "(*)i  (t^t)  n(.) 


CT  /■/-  -  1\  W1Z>1  _ \z  -  l)  _ 

‘  2  /  D(z>  (j-5-j-)  D(z)  +  KN(z) 
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with 


■p 

The  form  of  the  entries  in  the  column  vector  (W  GM) 1  (Eq.  61)  will 

1  9  * 

be  (ag7.  +  aiz  '  +  a2z  +  a3)/^(z)*  These  are  tabulated  in  Table  2  for  a 

basic  sampling  cycle  period  of  T  =  1.0  seconds.  Tne  first  column  corre¬ 
sponds  to  N(z)j  in  Eq .  61,  the  second  column  to  N(z)2»  and  so  on.  The 
feedback  gains  that  accomplish  all  design  objectives  are  tabulated  in 
Table  3,  The  design  objectives  are  repeated  below: 

•  Closed-loop  pole  assignment  as  specified  by 
Eq .  41. 

•»  Unity  steady-state  response  in  C  for  a  unit  step 
input  R. 

The  first  objective  is  satisfied  by  equating  like  coefficients  in  Eq .  41 
and  the  denominator  in  Eq.  64.  That  is, 

zD(z)  +  (z  -  l)KN(z)  =  A(z) 

=  (z  -  0. 1 353 )[ (z  -  0 . 1988)-  +  (0.3096)2] 

(65) 
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TABLE  2 


NUMERATOR  COEFFICIENTS  AS  A  FUNCTION  OF  ADVANCE  PARAMETER  A 


n 

0 

tt/10 

.4 

•  6 

.7 

a0 

3.000000000 

2.475580053 

2.251930038 

1.773697950 

1.584956596 

al 

-2. 3 1471762 1 

-1.623868245 

-1.340786259 

-0.756522207 

-0.537430945 

a2 

-0.195116512 

-0.029566508 

-0.095749095 

-0.206712756 

-0.235858883 

a  3 

-0. 067  316458 

-0.009062870 

-0.002312250 

-0.002620445 

-0.001415664 

Note:  D  (z)  =  (z  -  0. 367879441 ) (z2  +  0. 267796 1830z  +  0.018315639) 


TABLE  3 

.  FEEDBACK  GAINS 

K1 

0.091506404 

k2 

-1.494955697 

k3 

2.832023694 

k4 

-4.131755870 

K5 

2.703181447 
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0.4T,  0.6T,  and  0.7T  seconds.  These  individual  samples  are  stored  in 
appropriate  storage  registers  within  the  computer  until  all  samples 
within  a  group  are  obtained.  This  group  of  samples  is  then  added  toge¬ 
ther  with  each  sample  first  multiplied  by  an  individual  gain  factor  (Kj  , 
l<2 ,  K-j ,  ,  K^).  This  process  implies  that  a  full  frame  of  delay  equal 

to  T  seconds  is  inherently  present  in  every  component  except  the  first. 
That  is,  as  previously  pointed  out,  samples  cannot  be  used  to  form  the 
scalar  feedback  signal  in  Fig.  13  before  they  are  available.  This 
required  physical  delay  shows  up  in  the  analytical  model  as  Eq.  53 
(repeated  here): 


(KMdW]T  =  [KL,  K2z-1,  K3Z-1,  K4Z-’,  K^z"1]  (68) 

where  the  z~^  delay  factors  assure  that  samples  from  the  previous  frame 
time  (n  -  1)T  are  used  to  form  the  scalar  feedback  signal  for  the  n'f 
frame  time. 

H.  SOME  PROPERTIES  OF  THE  PSEUDO  MEASUREMENT  CONCEPT 

While  this  Is  not  a  particularly  realistic  design  problem,  a  number 
of  properties  of  the  pseudo  measurement  concept,  which  accrue  from  use 
ot  multiple-order  sampling,  have  been  demonstrated: 

•  Achieved  unity  steady-state  gain  without  forward 
loop  compensation. 

•  Generated  a  five  component  multiple-order  sample 
vector  from  a  scalar  analog  signal. 

•  Assigned  three  poles  using  only  the  one  available 
analog  measurement. 

•  Avoided  compensation  networks. 

•  Accomplished  pole  placement  with  multiple  samples 
of  a  single  output  variable  (in  distinction  to 
single  samples  of  several  output  variables). 

This  has  been  accomplished  without  using  an  ob¬ 
server  or  Kalman  filter  to  supply  estimates  of 
unmeasured  states. 
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We  remark  in  conclusion  that  the  pseudo  measurement  concept  is  not 
the  something-f or-nothing  gambit  that  it  might  appear  to  be.  In  connec¬ 
tion  with  the  illustrative  example,  a  price  is  paid.  This  price  is  more 
frequent  sampling  of  the  plant  output  variable  than  for  a  single-rate 
system  having  the  same  basic  sample  sequence  rate.  That  is,  additional 
samples  of  the  plant  output  variable  have  been  used  as  an  alternative  to 
using  digital  filters  (compensation)  or  as  an  alternative  to  measuring 
additional  plant  output  variables.  In  more  complex  multiple-order  sam¬ 
pling  schemes  the  general  nature  of  the  available  tradeoffs  will  be 
similar  to  the  alternatives  for  our  simple  example.  A  three-way  trade¬ 
off  (among  number  of  output  variables  measured,  number  of  samples  of  a 
given  set  of  measurements  taken  in  the  sampling  sequence  interval,  and 
complexity  of  digital  filtering)  will  exist.  Since  "number  of  output 
variables  measured"  translates  into  "number  of  sensors,"  the  tradeoff 
against  "number  of  samples"  greatly  favors  taking  additional  samples  in 
terms  of  practical  system  hardware  considerations. 

The  tradeoff  of  "number  of  samples  of  a  given  set  of  measurements 
taken  in  the  sampling  sequence  interval"  and  "complexity  of  digital  fil¬ 
tering"  also  favors  taking  additional  samples.  This  is  the  case  even 
though  the  number  of  registers  required  is  not  affected  strongly.  How¬ 
ever,  the  number  of  digital  multiplications  required  may  be  reduced  sig¬ 
nificantly.  This,  in  turn,  can  reduce  computational  burden,  which  is 
especially  Important  In  a  microprocessor  computing  environment.  Hence 
the  potential  contribution  of  the  pseudo  measurement  concept  to  cost 
reduction  in  system  hardware  is  realized. 

The  pse  ldo  measurement  technique  also  has  the  potential  to  provide 
for  gracefully  degrading  system  performance  in  the  face  of  sensor  fail¬ 
ure.  For  example,  if  a  sensor  failure  can  be  identified  and  isolated, 
then  it  is  possible  to  divert  the  computing  capacity  in  the  failed 
sensor  path  for  processing  of  additional  samples  of  different  output 
variables  having  operative  sensors.  For  this  use,  the  control  law  is 
restructured  in  a  predetermined  way  following  a  failure.  The  way 
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in  which  the  control  law  is  restructured  is  dependent  upon  the  particu¬ 
lar  failure.  This  may  enable  significant  reduction  in  number  of  redun¬ 
dant  sensor  components  when  coupled  with  so-called  analytical  failure 
identification  techniques.  The  latter  compare  the  measurements  of  dif¬ 
ferent  output  variables  for  consistency. 

Finally,  we  remark  that  the  sampling  sequence,  within  “he  basic  T 
frame  time  of  1,  e-T1//^s,  e~"^s,  e~'^s,  e~’,s,  was  selected  only  to 

demonstrate  a  multiple-order  sampling  environment.  There  is  nothing 
especially  magic  about  it;  in  fact,  it  might  even  be  a  poor  choice.  We 
cannot  tell  this  until  the  tools  needed  to  access  the  spectral  implica¬ 
tions  of  this  choice  (the  topic  of  the  next  two  sections)  are  developed. 


L: 
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SECTION  III 


FREQUENCY  RESPONSE  OF  SINGLE-RATE  SYSTEM 


A.  INTRODUCTION 

When  a  continuous,  stable,  stationary  linear  system  is  excited  by  a 
sine  wave,  its  steady-state  output  is  comprised  of  a  single  wave  at  the 
same  frequency  as  the  input-  It  differs  from  the  input  wave  only  by  a 
phase  angle  and  a  magnitude  factor.  Moreover,  it  is  unnecessary  to  com¬ 
pute  the  actual  transient  response  of  the  system  when  its  behavior  over 
large  epochs  of  time  is  of  interest,  since  both  the  magnitude  factor  and 
phase  angle  can  be  read  from  a  Bode  plot. 

A  similar  but  more  complex  situation  exists  when  a  sampler  and  data 
hold  couple  the  sine  wave  input  to  the  continuous  system.  This  configu¬ 
ration  i.3  referred  to  as  a  discretely  excited  cystem-  Civcn  that  the 
continuous  system  is  stable,  the  continuous  output  waveform  will  contain 
a  wave  at  the  fundamental  frequency  and  all  of  its  aliases.  Thus,  if 
the  system  is  forced  with  i  sin  bt ,  0  <  b  <  2w/T;  the  output  will  con¬ 
tain  terms  at  frequencies  b,  b  +  (2ir/T),  b  +  (4n/T),  ....  The  relative 
amplitudes  and  phase  angles  will  depend  on  the  data  hold  employed  as 
well  as  the  system  transfer  function.  Nevertheless,  given  the  data  hold 
and  system  transfer  functions,  the  magnitude  and  phase  angle  for  each 
and  every  component  can  be  read  trom  a  particular  "Bode  plot.”  Notice 
that  this  concept  of  frequency  response  is  more  comprehensive  than  the 
traditional  concept  of  the  "sampled  spectrum"  (see,  e.g.,  see  Ref.  3), 
which  is  limited  to  determining  the  single  sinusoid  that  fits  the  system 
output  samples  at  the  sampling  Instants. 

In  the  subsections  that  follow  we  review  the  frequency  response  con¬ 
cept  for  continuous  systems  and  then  proceed  to  the  frequency  response 
for  discretely  excited  open-loop  systems.  Finally,  the  manner  of  appli¬ 
cation  for  single-rate  sampled  closed-loop  systems  is  given. 
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B.  CONTINUOUS  SYSTEM  FREQUENCY  RF.SPONSF. 


It  will  be  helpful  to  first  review  the  frequency  response  eoneepL 
for  continuous  systems.  Let  R  in  Rig.  14  be  a  unit  sine  wave  input  with 
frequency  a  rad/sec.  The  Laplace  transform  of  the  output  Is 

C(s)  =  C(s)  R(s)  =  G(s)  (69) 

Sl  +  OJfi 


Figure  14.  Continuous  System 


Equation  69  can  be  expanded  in  partial  fractions  as 

gs  Terms  associated 

C(s)  =  — = — —  -y  +  -= — j+  with  characteristic  (70) 

S'  +  %  S'  +  ‘>)q  polynomial  of  G(s) 

The  bracketed  term  in  F.q.  70  determines  the  characteristic  behavior  of 
the  system  which  can  be  stable  (negative  eigenvalues)  or  unstable  (posi¬ 
tive  eigenvalues).  Nevertheless,  the  steady-state  sinusoidal  behavior 
is  completely  defined  by  the  partial  fraction  coefficients  A  and  13, 
since  once  they  are  known  the  steady-state  time  response  can  be  written 
directly  as: 

C(t)  .  ,  =  A  sin  M,,t  +  IS  cos  'J0t 

1  :1  m  u  u 

t+.« 

=  Va2  ■*-  h2  Bin  C-u0 1  +  v) 

where  ‘>  -  tan-'  (13/A).  The  details  of  solving,  for  A  and  1)  show  clearly 
the  relat loriship  between  the  Rode  plot  and  the  steady-state  waveform. 
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To  solve  for  A  and  B,  multiply  Eq.  70  by  [  s2  +  j  and  evaluate  the 
result  for  s  = 


G(s)“o|s-j»0 


or 


Terms  associated 
with  characteristic 
polynomial  of  G(s) 


(s2 


+  « 


21 
o i 


S“j“o 

(72) 


G(s) 


A  +  JB 


^A2  +  B2  eJ  tan-1  (B/A)  G(jw0) 


(73) 


To  summarize,  we  see  that  a  sinusoidal  input  at  frequency  wo  produces  a 
steady-state  sinusoidal  output  waveform  having  the  same  frequency.  It 
differs  from  the  input  only  by  a  magnitude  factor  and  a  phase  shift. 
Both  magnitude  factor  and  phase  shtft  for  any  given  input  frequency,  w0, 
can  be  read  directly  from  a  Bode  plot  for  G(jw0).  That  is,  for  any 
given  input  frequency  wQ, 


A  +  JB 


(74) 


The  sections  that  follow  expand  this  "frequency  response"  viewpoint  to 
include  dlsc-'etely  excited  systems. 


C.  MATHEMATICAL  PRKHMIHAR1KS 


Let  K  be  a  sinusoid  of  unit  amplitude  and  frequency  b  rad/scc.  If  K 
fs  sampled  at  1/T  samples  per  second  and  is  then  described  in  terms  of  a 
sample  sequence  having  N/T  samples  per  second,  the  result  is: 


k1' 


z^  yin  bT 

Z2N  _  2 (cos  bT)z^  +  l 


Z  “  c 


sT/N  * 


(75) 


*The  notation  follows  Ref.  1.  The,. superscript  denotes  the  period  oi 
the  sampling  operator,  f'ur  example,  K  and  K1'"1  denote  sampling  periods 
of  T  and  T/3  seconds,  respectively. 
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The  following  development  is  for  later  use.  First,  find  the  N  factors 
of  tlie  denominator  of  Kq.  75.  This,  in  turn,  permits  partial  fraction 
expansion.  This  expansion  contains  terms  which  correspond  to  sampled 
elementary  time  functions.  For  example,  if  f(t)  =  sin  bt,  then: 


i'T(t) 


[sin  bt]T  F(z )  =  — — - 


x  sin  hT 
2 (cos  bT)z  4- 


(76) 


However,  an  elementary  sampled  time  function  corresponding  to 


vr/.J  “  -9- - - -  (77) 

z2  +  2 (cos  bT)z  +  ) 

is  not  Immediately  apparent  (the  denominator  of  Kq .  77  is  one  of  the  N 
factors  of  the  denominator  of  Eq .  75).  This  situation  is  remedied  by 
addLof-  rt  to  the  trig  function  arguments  so  that  F(z)  in  Eq .  77  has  a 
recognizable  time  domain  counterpart.  lor  example,  Eq,  77  becomes: 


F(x) 


_ -x  sin  [  bF  (  h  /']')  ]T 

7?  -  2  cos  [ b+(*/T ) )T  z  +  1 


=£>  [f(t)JT 


-sin 


T 


(78) 


Consider  the  denominator  of  Eq.  75  for  the  special  case  of  N  =  2: 


2  cos  b'J’z2  4  1 


9/  bT\ 

-  /.I cos  —  Iz 


+  1 


4-  2 


(cos  H-) 


+  1 


(79) 


The  sign  of  the  coefficient  of  z  in  the  last  factor  on  the  right-hand 
side  of  Eq.  79  can  be  changed  by  adding  *  to  the  argument  bT/2: 


z4  -  2  cos  bTz2  4  1 


J 

/,  27T  \ 

T 

[COH 

( 

jb+T-) 

2 

4- 


(80) 
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Thus,  for  N  =  2,  the  partial  fraction  expansion  contains  terms  which  corre¬ 
spond  to  sampled  sine  and  cosine  functions  not  only  at  the  input  frequency, 
b,  but  also  at  the  first  alias,  b  +  (2ti/T),  of  the  input  frequency. 

f 

The  general  result  is  obtained  by  induction  by  repeating  the  exercise 
for  N  =  Nj  and  for  N  =  Nj  +  1  (Ref.  l)i 


(82) 


Thus,  the  sampled  sine  wave,  [sin  btjT,  can  be  described  alterna¬ 
tively  in  terms  of  a  sum  of  sinusoidal  components  at  frequency  b  rad/sec 
and  its  first  N  -  1  aliases  sampled  with  period  T/N. 

D,  OPEN-FjOOP  FREQUENCY  RESPONSE  — -  CONTINUOUS  OUTPUT  SAMPLED 

Consider  the  system  of  Fig.  15  where  G(s)  represents  an  arbitrary 
transfer  function  and  M  represents  an  arbitrary  data  hold.  Suppose  R  is 
a  unit  amplitude  sLne  wave  and  the  output  is  sampled  with  period  T/N. 
Using  multi-rate  sampling  results  from  Section  II, 

C  '  /»  =  (  gmpt ‘ l  /  N  -  <CM)T/»  kt 


-  <CM)T/«  ~.2.- 


sin  bT  _ 

2 (cos  bT)zN  +  1 


z  =  cBT/N 


(83) 
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Expand  L he  right-hand  aide  of  Kq .  83  in  partial  fractions: 


C'f/M 


A,-,z  sin  un(T/N)  +  Knz[z  -  cos  (^(T/N)] 
z^  -  27.  cos  «n(T/N)  +  1 


+  (Terms  due  to  modes  of  (GM)T/hT]  (84) 


R  /  R* 

M 

C  y  CT/N 

T 

T/N 

Figure  15. 

Open 

-loop 

System 

Assume  that  responses  in  the  modes  of  (GM)^N  approach  zero  as  t  +  », 
i.e.,  that  all  modes  are  stable.  This  restriction  simplifies  the  nota¬ 
tion  in  the  derivation  (the  development  that  follows  Is  valid  whether 
the  modes  of  (CM)  are  stable  or  unstable).  In  Eq .  84, 

w„  =  b  ,  n  =  0, 1,2,.. .  ,N-1 


For  the  present  ,  assume  b  <  2n/T.  The  steady-state  waveform,  at  the 
sampling  instants,  can  be  written  as: 


;  c(  t  >  ] 


T/N 


N-l  V/N 

2  (An  sln  wnt  +  Bri  cos  wnL) 
n=0 


(85) 


To  solve  for  A  and  li^ ,  multiply  each  side  of  Eq .  84  by 

!  z -  2z  cos  «  (T/N)  +  l]  ,  0  <  k  <  (N-l) 

and  evaluate  ‘or  v.  =  (T/N) .  The  only  non-zero  term  on  the  right- 

hand  side  of  Lh<'  result  occurs  for  n  =  k . 


(gh)T/N  -2"  «in_y  , 


jz — z~nr\7-  -2  003  a-k(,J'/N)z  + 1 ) 

( ?  cos  d*i )  vr  +1 


1/^('L'/N) 


[Anz  sin  eun(T/N)  1  +  Bnz[z-oos  cy/T/ft)  ]  0 

2->  ■  ;  7~7, - - -  lz‘  ~Sz  cos  ac(T/N)  -•  1  ] 

nO  E  cos  o^Cr/Hjz  +  1  * 


For  any  n  f-  k,  the  right-hand  side  of  Eq.  86  is  identically  zero  since 
z2  -  2[cos  a)  (T/N)]z  +  1  =  [z  -  cos  to  (T/N)]2  •(-  (sin  w  (T/N)]2 


vanishes  when 

z  =  Uwr(T/N)  -  cos  ur(T/N)  +  j  sin  wr(T/N)  (88) 

Evaluation  of  Eq.  86,  subject  to  the  constraint  of  Eq.  88,  is  a 
tedious  chore  (Ref.  1).  However,  the  result  is  relatively  simple: 


+  JB"  ■  T,  <™>T/\,14„n 


=  psT/n 


(T/N) 


Notice  in  particular  the  definition  of  z  used  in  the  evaluation. 

To  summarize,  when  the  system  is  forced  by  sin  bt ,  the  steady-state 
output  waveform,  sampled  with  period  T/N,  is  given  by  Eq .  85.  The  coef¬ 
ficients  An  and  Bn  in  Eq .  85  are  computed  using  Eq .  89. 

For  example,  let 


1  -  e~sT 


G(s)  = 


s  +  a 


so  that 


(91) 


~  <GM)T/N 


1  -  e-aT/N  (1  -  z~N) 
N(z  -  e“aT/N)  (l  -  z~l) 


It  is  instructive  to  plot  the  amplitude  ratio  of  the  frequency 
response  for  Eq.  91  with  N  as  a  parameter.  For  the  sake  of  clarity,  the 
plot  is  versus  w  rather  than  log  oj. 

Refer  to  Fig.  16.  Notice  that  over  the  range  of  frequency  for  the 
plot,  0  <;  a)  <  8ti/T,  the  N  =  1  case  amplitude  ratio  folds  about  the  fre¬ 
quencies  tt/T,  3n/T .  7tt/T;  the  N  =  2  case  folds  about  frequencies 

2n/T,  4tt/T,  6ti/T;  and  the  N  =  4  case  folds  about  m  =  4ir/T.  However,  it 
can  also  be  observed  that  the  various  amplitude  ratio  plots  are  periodic 
with  frequency.  For  N  =  1  this  period  is  2ir/T;  for  N  =  2,  4n/T  and  for 
N  =  4,  8rr/T.  Each  period  cotitains  precisely  the  number  of  frequency 
points  (at  the  input  or  its  positive  alias  frequencies)  required  to 
match  the  continuous  steady-state  time  response  at  the  sampling  instants 


Figure  16.  Magnitude  Plot  for  N  =  1,  2,  4 
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and  at  N  -  1  equally  spaced  inter-sample  points.  Consequently,  It  Is 
the  sampling  frequency  Interval  rather  than  the  folding  frequency  inter¬ 
val  that  is  fundamental  to  generalizing  the  concept  of  frequency  re¬ 
sponse  . 

Consider  the  use  of  Fig.  16.  Imagine  that  a  unit  amplitude  sine 
wave  with  a  frequency  it / 2  is  input.  In  the  N  =  1  case,  our  interest  is 
matching  the  continuous  steady-state  time  response  only  at  the  sampling 
instants  with  a  single  sine  wave.  The  magnitude  and  phase  angle  (not 
shown  in  Fig.  16)  can  be  read  from  this  plot  at  to  =  n/2  (or  n/2  +  2tt/T, 
n/2  +  4nT,  rr/2  +  6ir/T,  ...;  any  of  these  points  gives  the  correct 
valu^).  Clearly,  if  the  objective  is  to  match  at  the  sampling  instants 
with  a  single  sinusoid,  the  frequency  could  be  b  plus  any  2n/T  multi¬ 
ple.  One  cannot  tell  the  difference  once  the  waveform  is  sampled.  In 
fact,  the  "sub"  aliases  at  b  -  (2n/T),  b  -  4n/T  will  also  work.  These 
"suh"  aliases  are  the  difference  terms  which  are  prominent  in  modular Luu 
theory. 

The  N  =  1  plot  in  Fig.  16  corresponds  to  the  "sampled  spectrum"  fre¬ 
quency  response  of  sampled  data  control  theory.  Consider  the  N  =  2  case 
wherein  the  objective  is  to  match  one  inter-sainple  point  as  well  as  the 
sample  points.  Let  the  input  freuency  be  n/2  and  note  that  the  points 
at  u  =  n/2,  n/2  +  2n/T  give  the  magnitude  (and  phase)  for  the  sine  waves 
at  those  frequencies,  as  would  the  points  n/2  +  4rr / T ,  tt/2  +  6ti/T.  Sup¬ 
pose  next  that  the  input  frequency  is  b  =  n/2  +  2n/T.  Clearly,  the 
second  required  component  could  be  read  from  the  "first  alias"  at 
b  +  2tt/T  or  the  first  sub-alias  at  b  -  2n/T  (or,  for  that  matter,  at  a 
host  of  other  frequencies). 

In  the  N  =  4  case,  four  sine  waves  are  required  to  fit  three  inter¬ 
sample  points  as  well  as  the  sample  points.  If  the  Input  frequency  were 
b  =  n/2  +  6n/T,  and  if  the  plot  of  Fig.  16  with  its  limited  range  of  8n 
were  the  only  one  available,  clearly  it  would  be  to  our  advantage  to  use 
the  "difference"  frequency  points  at  to  =>  b  -  2n/T,  b  -  4n/T,  and 
b  -  6n/T  to  establish  the  magnitude  and  relative  phase  of  the  other 
three  sine  waves. 
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Consider  now  the  generalization  of  these  observations.  Notice  that 
the  N  s l lie  wave  components  required  to  match  the  continuous  steady-state 
time  repor.se  at  the  sampling  instants  and  at  N  -  1  equally  spaced  inter¬ 
sample  points  need  not  have  frequencies  from  within  one  period  of  the 
amplitude  ratio  plot.  However,  the  frequencies  used  must  each  be  an 
alLas  separated  by  (2ir/T)Nk  (where  k  is  an  arbitrary  integer  for  each 
frequency)  from  the  N  alias-related  frequencies  falling  in  any  one 
period  of  the  frequency  response.  This  is  evident  in  the  example  for 
N  =  2.  If  the  input  frequency  Is  tt/2,  one  would  ordinarily  choose  sine 
wave  components  with  frequencies  tt/2  and  [tt/2  +  (2if/T)J  to  represent 
the  continuous  steady-state  response  at  the  sampling  instants  and  one 
inter-sample  point.  Alternative  choices  of  frequencies  are  tt/2  and 
[tt/2  +  (6it/T)J,  and  so  on.  Choices  of  n /2  and  [tt/2  +  (4tt/T)]  or 
[tt/2  +  (2ir/T)[  and  [tt/2  +  (6ti/T)[  are  not  alternatives  because  these 
frequency  choices  do  not  satisfy  the  separation  criterion  stated  above. 

This  brief  discussion  associates  aliases  and  sub-aliases  with  the 
sum  and  difference  frequencies  of  modulation  theory.  It  is  not  the 
case,  however,  that  sum  and  difference  components  are  necessarily  both 
present  simultaneously  in  the  output.  It  is  shown  above  that  only 
N  components  are  needed. 

The  case  where  N  is  extremely  large  Is  also  of  interest.  Let  N  +  «> 
after  evaluating  Eq.  91  at  z  =  1  u)n(T/N): 


1  _  e-aT/N 
n(z  -  e_aT/NJ 


i  -  z_N! 

i—T-Tl 

lz-14u>„(T/N) 


_ [l  -  e~aT/N)(l  -  1  -ai nT) _ 

N[l  wn(T/N)  -  e-aT^NJ[l  -  1  -«n(T/N)] 


lira 

N+w 


(92) 
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An  Indeterminate  form  is  obtained.  Therefore,  use  L'Hopital's  rule 
twice  [substitute  1  <sn(T/:})  -  cos  wn(T/N)  +  j  sin  wn(I/N),  etc.]  and 
obtain  (see  Appendix  I  of  Ref.  1): 


1  1  -  e 


-sT 


N  s(s  +  1) 


z  =  Uwn(T/N) 

1  -  e'j01«T  1 


lim 
]  N+® 


j“nT  1  +  3wn 


1  -  e~sT 

sT  s  +  1 


(93) 


s.i“>n 


where  z  =  esr,/N.  That  is,  as  N  ■>  “>,  simply  divide  GM  by  T  and  evaluate 
the  coefficients  at  s  -  juin.  This  is  representative  of  the  general 
result  of  the  next  subsection. 

E.  OPEN-LOOP  FREQUENCY  RESPONSE:  CONTINUOUS  OUTPUT 


In  the  previous  subsection  it  was  shown  that 


*n  +  3®n  "  1 


(T/N) 


£  P3T/N 


z  =  e 


(94) 


To  deduce  the  behavior  for  infinite  N,  rewrite  Eq.  94  as 


f i!  -  I  Jl  <«,(■♦  3^)]} 


(95) 


The  T/N's  cancel,  and  as  N  gets  very  large  only  the  k  =  0  erm  contri¬ 
butes  since  all  "aliased"  components  of  the  spectrum  have  moved  to  in¬ 
finite  frequency.  Therefore, 


44 
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i™  5  <™>T/NU-L  „„(!/») 


Z 


(96) 


GM 

T 


|s  =  jun 


esT/N 


2  un 


wn  =  b  +  T 


where 


l~\ 

V'WlNT 


n  =  nG ,  riQ  +  1 ,  . . .  ,  N  -  n0  -  1 


2tt 

T 


The  example  of  the  previous  section  with  N  finite  can  now  be  studied 
for  N  infinite.  This  result  gives  the  frequency  response  for  the  con¬ 
tinuous  output. 


An  +  jBn 


1 


s  +  1 


s  =  jwn 


(97) 


A  Bode  plot  for  this  result  is  shown  in  Fig.  17.  Components  for  the 
input  frequency  b  =  1.0  rad/sec  and  its  aliases  have  been  indicated  with 
the  square  symbols. 

Interpretation  of  Fig.  17  is  as  follows.  Suppose  a  unit  sine  wave 
at  1  rad/sec  is  input  to  the  sampler.  Then,  if  the  amplitude  ratios  and 
phase  angles  for  sine  waves  with  frequencies  1,  (1  +  2tt/T),  (1  +  4m /T), 
....  are  read  from  Fig.  17,  and  these  components  are  added  together,  the 
resultant  waveform  will  be  an  exact  natch  of  the  actual  steady-state 
output  waveform  in  Fig.  18.  One  might  expect  this  waveform  to  be  very 
nearly  sinusoidal,  since  the  first  alias  is  attenuated  on  the  order  of 
30  dB,  relative  to  the  input  component.  However,  the  steady-state  time 
response  does  not  bear  this  out,  as  can  be  seen  in  Fig.  18.  The  reason 
is  that  the  summation  of  alias  terras  is  significant  despite  their  small 
Individual  size. 

Again,  in  connection  with  Figs.  17  and  13,  consider  a  (1  +  4m/T) 
rad/sec  unit  sine  wave  input  to  the  sampler.  The  steady-state  output 
will  contain  sub-aliases  at  1  rad/sec  and  (1  +  Zm/T)  rad/sec;  a  compo¬ 
nent  at  the  input  frequency,  (1  +  4m /T)  rad/sec;  and  other  aliases  at 
(1  +  6m /T),  (IT  8m /T),  rad/sec.  The  amplitude  ratio  and  phase 

angle  for  each  component  is  again  read  from  Fig.  17  (at  points  indicated 
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bv  tlie  square  symbols).  VJhen  these  components  are  added  the  steady- 
state  time  response  is  again  that  in  Fig.  18.  Thus,  the  continuous 
response  spectrum  and  steady-state  time  response  are  the  same  for  a  sine 
wave  input  of  given  frequency  as  for  a  sine  wave  input  having  a  positive 
frequency  which  is  an  alias  of  the  given  frequency. 

Another  observation  is  that  the  aliases  do  not  represent  "harmonic" 
components  but  rather  represent  modulation  components  which  must  add 
together  properly  in  order  to  match  conditions  at  the  sampling  instants 
for  the  haste  period,  T.  It  can  be  seen  in  Fig.  18  that  the  "steady 
state"  does  not  necessarily  imply  "periodic."  Periodic  waveforms  occur 
only  when  the  input  frequency  and  the  sampling  frequency  have  an  integer 
relationship  with  respect  to  one  another. 


F.  SINGLE— RATE  CLOSED- LOOP  FREQUENCY  RESPONSE 

Closed-loop  results  are  dependent  upon  the  conf igurat ion  of  the 
digital  loops.  However,  the  basic  analysis  procedure  Is  independent  of 
configuration.  It  is  important  to  understand  this  procedure,  and  in 
particular  the  simplifications  that  occur  in  analysis  of  closed-loop 
systems . 

Consider  the  (vector)  system  shown  in  Fig.  19.  The  objective  is  to 
find  the  coefficients  that  characterize  the  spectral  components  of  the 
continuous  response,  C,  i.e.,  the  Frequency  response.  The  proce<u  re  for 
this  example  is  typical.  First,  solve  for  the  vector  component  at  the 
input  of  tiie  data  holds. 
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fJ  =  -  GTGT(GM)teT  (98) 

Theref  ore , 

ET  =  [I  +  GTGT(GM)T]_1  GTrT  (99) 

Solve  for  C(s)i 

C  -  (CM)  [I  +  GTGT(GM)t]"1  GTrT  (100) 

The  spectrum  of  C(s)  Is  of  interest.  It  is  determined  by  first  finding 
the  spectrum  of  and  taking  the  limiting  case  as  N  + 

Let  the  input  be  a  unit  sine  wave  at  frequency  b  rad/sec  and  let 


,sT/N 


(101) 


so  that 


RT  = 


sin  bT 


Z2N  _  2zL  (cos  bT)  +  1 
Using  Eq.  102  and  the  results  of  Se.ction  II, 


(102) 


cT/n  =  (gm)t/n[i  +  GTg|(gm)t]-1  g|rt 


(103) 
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For  the  sake  of  brevity  write  Kq .  103  as 


c T  /  f-1  =  GT/0  gTrT 

A  R 


(104) 


Expand  the  right-hand  side  of  Eq .  104  in  partial  fractions: 


rT/N  -  rT/M  rT _ 7.  1  sin  bT _ 

A  B  Z2N  _  2z N  cos  p,-p  4.  i 


y]  Anx  sin  un(T/N)  +  B11z ( z  -  cos  uin(T/N)] 
n  —  q  -  2z  cos  <jJn(T/N)  +  1 


+  [Terms  due  to  modes  of  gJ/^  G^l 

A  D 


(105) 


Assume  that  responses  in  the  modes  of  g]/N  q!  approach  zero  as  t  M. 
t.c.,  that  all  closed-loop  system  modes  are  stable.  In  Eq.  105, 


wn  =  b  + 


n  =  0,1,2,..., N-l 


(106) 


Notice  that  Kq .  105  is  exactly  the  same  as  Eq .  83,  except  (GM)T,/N 
lias  been  replaced  by  G J  Hence  a  crucial  result  is  obtained  using 

Eq.  89: 


An  +  pn  “  b  GI/NgI1 

r  N  A  B|z=Uwn(T/N) 


(107) 


Cg(z )  =  CT(ZN) 


z  =  es^ 


008) 


because  of  the  local  definition  of  z  as  es^/^.  Tiierefore,  using 
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I 


1  Uwn(T/N]N 


I 

U«nr 

cos  (i)nT  +  j  sin  fOnT 
cos  (b  +  (2*n/T) ]T 

+  j  sin  [b  +  (2'irn/T)  ]T 
cos  bT  +  j  sin  bT 


we  obtain 


Uwn(T/N) 


Gl(2) 


?.=!4bT 


(109) 


Thus,  z1'1  in  can  be  replaced  with  z  provided  it.  is  evaluated  for  z  ~ 
1 4bT  instead  of  for  z  -  14<d  (T/N).  Hits  amounts  to  a  simple  (and  very 
convenient)  change  of  variable.  At  this  point  the  result  is: 


An  + 


jRn  “  4  Gj/N(z)  g'|(/;) 

11  N  A  ?.-liu)n(T/N)  B 


z- 1  Jf.  bT 


=  esT/N 


(HO) 


Equation  110  is  the  basic  result  for  the  finite  N  case.  To  reiterate, 
to  find  the  coefficients  of  the  N  sine  waves  matching  the  T/N  sampled 
output,  C,  compute  the  usual  pulsed  transfer  functions  for 

GjCz)  -  f  I  +  c[G|(GM)TrlGj  v  z  =  esT  (111) 
and  evaluate  it  for 


z  =  l$bT  (112) 

Next,  compute  the  usual  T/N  pulsed  transfer  function  for  as  a  func¬ 
tion  of  z  and  evaluate  it  for  z  =  ltfw^O'/N)  where.  wn  •=  b  +  (2irri/T). 
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Thus,  it  is  superfluous  to  consider  g£  as  a  function  of  z^;  it  suf- 

o 

fices  to  consider  it  as  a  function  of  2.  Moreover,  only  G^/N  is  a  func¬ 
tion  of  N;  this  simplifies  obtaining  the  limiting  case  tremendously. 
For  the  case  of  N  +  the  continuous  case,  we  obtain 


An  +  J^n 


M(s)GA(s) 


s-jmn 


)( 


GT(z) 


,sT 


z  =  UbT- 


(113) 


Equation  113  is  the  desired  result  for  the  given  closed-loop  config- 
uration.  However,  it  is  the  procedure  which  is  key.  One  can  follow 
the  details  through  for  other  configurations  quite  easily. 

An  illustrative  closed-loop  frequency  response  example  is  shown  in 
Fig.  20. 

G.  SUCTION  SUMMARY 

The  "sampled  spectrum"  frequency  concept  of  sampled  data  control 
theory  is  concerned  with  determining  the.  single  sinusoid  that  fits  the 
output  of  a  single-rate  system  at  the  sampling  instants.  In  this  report 
the  frequency  response  concept  has  been  extended  to  encompass  the  con¬ 
tinuous  spectrum  for  the  continuous  variables  of  a  discretely  controlled 
system.  Moreover,  the  theory  is  sufficiently  comprehensive  to  cover 
cases  wherein  a  group  of  N  sinusoids  is  used  to  roatch  the  continuous 
variables  not  only  at  the  sample  points  but  at  N  -  1  equally  spaced 
inter-sample  points  as  well.  N  may  be  finite  or  infinite.  Infinite  N 
corresponds  to  the  true  continuous  spectrum.  The  practical  value  of 
knowing  the  true  coutinuous  spectrum  for  continuous  physical  variables 


*Ac.curate  numerical  determination  of  Gj^(z)  |  z = j  ^  i^T  raav  Prove  diffi¬ 
cult  at  high  sampling  rates.  This  is  the  result  of  small  differences 
between  large  numbers  which  occur  in  the  computations  as  poles  and  zeros 
approach  the  unit  circle.  In  this  event,  one  is  well  advised  to  carry 
out  equivalent  computations  in  a  domain  where  numerical  conditioning  Is 
much  Improved  (e.g.,  in  terms  of  w'  or  w). 
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Figure  20.  Magnitude  Plot  for  Closed-Loop  bystem 


is  self-evident.  For  example,  this  information  can  be  used  for  evaluat¬ 
ing  inter -sample  ripple  and  for  guiding  continuous  filter  selection. 
Cases  wherein  N  is  finite  also  have  important  practical  application  as, 
for  example,  in  the  bench  validation  of  the  digital  controller  hardware/ 
software  combination.  Digital  controller  characteristics  are  often 
specified  in  terms  of  end-to-end  "frequency  response."  On  the  bench,  if 
continuous  outputs  (say  from  data  holds)  are  sampled  at  a  finite  rate  (N 
is  finite),  the  difference  between  the  measured  frequency  response  and 
the  true  continuous  frequency  response  (N  +  “>)  may  be  significant.  The 
results  of  this  report  can  be  used  to  predict,  minimise,  or  correct  for 
the  difference  between  measured  and  true  frequency  responses. 

, Results  for  closed-loop  cases  depend  upon  the  specific  configuration 
of  the  digital  loops.  However,  the  basic  analysis  procedure  is  indepen¬ 
dent  of  the  loop  configuration.  The  procedure  is  to  obtain  the  pulsed 
transfer  function  matrix  relating  the  outputs  to  the  data  holds  and  the 
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inputs  through  which  the  (sampled)  sine  wave  inputs  enter.  The  continu¬ 
ous  system  (including  data  hotds)  is  placed  in  cascade  with  the  output 
of  the  pulsed  system.  Frequency  response  is  then  evaluated  in  exactly 
the  same  manner  as  for  any  discretely  excited  open-loop  continuous 
system. 

Corresponding  results  for  two  rate  sampled  closed-loop  systems  are 
contained  in  Ref.  1.  Extensions  for  true  multi-rate  sampled  closed-loop 
systems  will  be  developed  in  Section  IV.  In  addition,  the  theory,  in 
Section  IV,  is  extended  to  provide  a  tool  for  quantifying  fidelity  of 
digital  simulations  of  continuous  systems. 
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SECTION  IV 


FREQUENCY  RESPONSE  OF  MULTI-RATE  SYSTEMS 


A.  INTRODUCTION 

The  frequency  response  of  digitally  controlled  single-rate  systems 
was  discussed  in  the  previous  section.  The  manner  in  which  one  treats 
two- rate  systems  is  given  in  Ref.  1.  We  will  omit  a  review  of  the  two- 
rate  case  and  move  directly  to  a  discussion  of  multi-rate/multiple-order 
configurations,  it  will  be  shown  that  the  basic  tools  developed  in  this 
section  are  ideally  suited  to  a  variety  of  situations  such  as  defining  a 
Bode  plot  for  digitally  controlled  closed-loop  systems,  a  Bode  plot  for 
the  simulation  of  a  closed-loop  system,  and  the  Bode  plot  associated 
with  a  simulation  wherein  more  than  one  independent  processor  (e.g., 
multiple  computers  or  microprocessors)  is  used  to  implement  the  computer 
code. 

B.  A  BASIC  RESULT 

Let  the  general  nrul ti-ra te/multiple-order  open-loop  system  of  Fig. 
21  have  a  sine  wave  input. 

Tn  Flo.  21 

- - iP'  - f 


CB  *  ( G i Ra ) ^  (114) 

CT/N  =  [G[GlRa)0]T/N  (115) 

where  a,  6  represent  sampling  schemes  with  a  basic  period  of  T 
seconds.  For  example,  a  might  represent  a  multiple-order  sampling 
format;  B  might  represent  a  multi-rate  and/or  pseudo  measurement  format 
(Section  II). 


Figure  22.  Open-Loop  System  with  Switch  Decomposition 


If  (i  represents  multiple-order  sampling  ana  n  a  pseudo  measurement 
format  usinp  multi-rate  sampling,  the  switch  decomposition  modeling, 
components  might  appear  as: 


5^ 


(118) 


Wj  =  fl,  e_T°S,  e'TlS,  e~T2S; 


'1* 


1 

Jo3 

2Tis 

qt2s 


(119) 


and 


W  2  =>  W  2  = 


0  0 

e-sT/3  0 

0  e'2sT/3 


(120) 


J2* 


1 

rn  /  o 

gSii  /  J 

P2aT/3 


(121) 


Comparing  Eq.  117  with  Eqs.  83  (or  Eq.  104),  we  see  that  only  one  new 
facet  has  entered  the  problem,  namely  (W,  R)T  replaces  RT.  Consider, 

rp 

therefore,  a  generic  component  of  (Wj  R)  -  for  instance  the  scalar 

(e^sR)^,  where  0  <  6  <  1.0. 

For  R  =  sin  bt  and  A  zero  (e^s  -  1),  the  output  equation  (Eq.  117) 


becomes 


CT/N  =  (GW2)T/N(W2*G1W1)T 


sin  bT 


Z2N  _  (2  cos  bT)z^  +  1 


=  pST/N 


(122) 


T 

where  Rl  is  described  in  terms  of  a  N/T  samples  per  second  model.  For 
the  sake  of  brevity  write  Eq.  122  as 


CT/N  _  CT/N  cT  _ zN  sin  bT1 _ 

A  B  Z2N  _  (2 

cos  bT)zll  +  1 


(123) 


56 


Expand  the  right-hand  side  of  Eq.  123  in  partial  fractions: 


(-T/N  =  gT/NGT  zN  sin  bT 

A  B  Z2N  _  (2  cos  bT)z^  +  1 

Anz  sin  wn(T/N)  +  Bnz[z  -  cos  un(T/N)] 
z2  ~  [2  cos  oin(T/N)]z  +  1 

+  [Terms  due  to  modes  of  (124) 

For  non-zero  A,  we  use  the  advanced  z-transform  on  (e^sR)^  and 
Eq.  124  becomes 


CT/N 


gT/NgT  zN[(sin  bAT)zN  +  sin  [b(l  -  A)T ] 
A  B  z2N  _  (2  cos  bT)zN  +  1 


N-  1 

r 

n=0 


Anz  sin  u)n(T/N)  +  B^tz  -  cos  wn(T/N)] 
z2  -  [2  cos  un(T/N)]z  +  1 


+  [Terras  due  to  modes  of  G^^Gg]  (125) 

Assume  that  responses  in  the  modes  of  G^^G^  approach  zero  as  t 
i.e.,  that  all  modes  are  stable.  In  Eqs.  124  and  125 

(un  =  b  +  —•  ,  n  =  0,1,2 . N-l  (126) 

(More  generally,  n  =  n..,  n_+I. ,  N-n  -1,  where  n_,  =  -(b/a>„)  ].  The 

o  o  o  os  INT 

steady-state  waveform,  at  the  sampling  Instants,  can  be  written  as 


C(t) 


T/N 


N-l 


2  (An  sin  wnt  +  Bn  cos  wnt) 
n=0 


T/N 


(127) 
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To  solve  for  ^  and  Bn,  multiply  each  side  of  Eqs.  124  and  125  by 

{z2  -  [2  cos  uk(T/N)]z  +1}  ,  0  <  k  <  (N-l)  (128) 

and  evaluate  for  z  =  1  «k(T/N).  Since  the  only  terms  which  can  survive 
on  the  right-hand  side  occur  when  n  =  k,  the  £  will,  of  course,  dis¬ 
appear.  Then,  if  we  so  choose,  the  k  notation  can  be  changed  back  to  n. 
To  illustrate. 


“A  “B  2N  - |  *2-  U  cos  w  (T/N)  ]z  +  1 } 

A  B  z  ■  <2  cos  bT)z%  1  '  >z=Mu)k(T/N) 


gT/ngt 


V*1  Anz  sin  Wn(T/N)  +  Bnz[z-a.s  «>„(T/N)]  |?  ,,  ,  J 

Y - - z'-  [2  cos  iok(T/N)]z  +  ] . 

n*Q  .  —  m  /lt  /t/mi  l. ,  i  '  * 


z  -(2  cos  wn(T/N)lz+l 


z«Mi.fc(T/N) 


(129) 


qT/NqT  «,?  +  ,«  in  b '  1  -  fl)T  li  z2  _  ;  2  CpS  M  (T/N))z+  ij 

A  B  z2N-(2  cos  bT)zN  +  l  1  K  1 


z-12tak(l'/N) 


N-l 

E 

n'O 


Apz  sin  u^,(T/N)  +  Bnz(z  -  cos  wn(T/N) )  |  ,  i 

Y.  — - - - — - 2 -  z2-12  cos  wk(T/N)]z+l 

z  -[2  cos  wn(T/N)]z+l  1  ’ 


(130) 


z=14wk(T/N) 


For  any  n  j4  k,  the  right-hand  side  of  Eqs.  129  and  130  are  identically 
zero  since 


z2  -  [2  cos  o>k(T/N)]z  +  1  =  [z  -  cos  u>k(T/N)]2  +  [sin  «k(T/N)  ]  2 

(131) 


vanishes  when 

z  =  ]£wk(T/N)  =  cos  «k(T/N)  +  j  sin  ^(T/N)  (132) 
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Specifically,  we  obtain 


[cos  w  (T/N)  +  .1  sin  <u  (T/ll)  -  cos  ^(T/N)]2  +  [sin  <i'k(T/N)]2  =  0 

(133) 

For  n  =  k,  the  cancellation  of  the  common  factor  guarantees  the  survival 
of  an  n  =  k.  term.  Factoring  out  a  common  z  gives 

z(Ak  sin  w  (T/N)  +  B^fcos  w  (T/N)  -  cos  w^CT/N)  +  j  sin  w^CT/N)]} 

=  (A.  +  jB.  )z  sin  0).  (T/N)  I 

^  J  k  k  |z«l*u  <T/N) 

(134) 

Tlicrefore,  Eqs.  129  and  130  become 


/■  sill  u^d/H)  +  jbk) 


cI,NcI(^ 


sin  bT)jz2-  [’cos  mk(T/H)]z  +  1  j 
z2N-  ( 2  cos  bX)zN  +  1 


z-Uuik(T/N) 


(135) 


sin  iu  (T/N)(Ak  +  J»k> 


?/  |vsln  i>AT)  +  sin  b(  1  —  ii)  T  ]  j  z2  -  [7  c.os  0Jk(T/N)]  z-l-  1  j 


cos  bT) +  1 


z-i^ix/n) 


(136) 


At  this  point,  divide  through  by  z  sin  u)^(T/N)  and  replace  k  with  n: 


A„1  JB„ 


A  1! 


c 


z-l-.U)n(T/N) 


zN~1  sin  bT 
sin  0^,(1  /U) 


z2-  1 2  cos  Wn(T/N)J  r.  H  l| 
z  -  (2  cos  l>T)  zN  +  1  j 


037) 


An  ■*  J  I'n 


,T/Nr'l 
k'A 


*«L'.u^(T/N) 


zN-’  [(sin  bAT)zNl  sin  b')-6)Tl 
sin  U^CT/N) 


z?-  -  [2  ian(T/N)]z  +  1 

z ^  —  (2  cos  bT)zN4  1 


(T/N) 

n 


(138) 
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The  last  term  on  the  right-hand  side  of  Eqs.  137  and  138  is  indetermin¬ 
ate  (0/0)  when  z  =  ljfoJn(T/N).  Therefore,  applying  L'Hopital's  rule 
first  to  Eq.  137  (A  -  0  case), 


An  +  JBn 


GT/NGT 
A  B 


zN  1  sin  bT  2z  -  2  cos  ni^CT/N) 
sin  wn(T/N)  2N[z^N~1  _  (COs  bT)z^-^) 


z=l<«n(T/N) 


(139) 


GT/NGT  sin  bT  2[z  cos  h)n(T/N)] 

A  B  sin  o'n(T/N)  2Mz^~^[zN  -  cos  bT)  z=i^w  (T/N) 


(140) 


An  +  JBn 


gt/ngt 

A  B  z=l  uin(T/N) 


sin  bT[cos  wn(T/N)  +  j  sin  ton(T/N)  -  cos  wn(T/N)] 

[ sin  un(T/N)  ]N  j  [  Ucon(T/N)TN  "  cos  bT)}  ^ 


An  +  JBn 


°A7 NC b) z=14wn( T/N) 


(sin  0T)[j  sin  h>n(T/N)] _ 

[sin  un(T/N)]N j[(l4wn(T/N)]N  -  cos  bT J 


but. 


[Uwn(T/N)JN  -  cos  bT  =  l*wnT  -  cos  bT 


(142) 


=  cos  wnT  +  j  sin  wnT  --  cot  bT 

=  cos  [b+(2irn/T)  ]T  +  j  sin  [b+(2iTn/T)]T'  -  cos  bT 
=  cos  bT  +  j  sin  bT  -  cos  bT 
=  j  sin  bT 
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Therefore, 


An  +  JBn 


I  gt/ngt 

N  A  B 


z=Uun(T/N) 


(143) 


Next , 


G}<*)  =  c|(zN) 


z  =  eST 


(344) 


because  of  the  local  definition  z  =  esT/N.  Therefore, 


GT(ZN)|  =  GT[U«n(T/N)]N 

V  |z=Uu.n(T/N)  B 


=  G'J(1*u„T)  =  Gj(UbT)  (1«) 


l-hat  is?,  take  the  "T"  z-transform  of  G?,  and  evaluate  at  z  -  WbT. 
The  finite  N  result  for  Eq.  137  (A  =  0)  is  then 


An  +  jBn 


gT/n 

A 


z  ^esT/N 

GR 

z=U<an(T/N)  B 


a 

z=es  L 
z=UbT 


(146) 


As  M  +  «,  the  coefficients  of  the  continuous  spectrum  are  given  by 


1  I  ~Tlz-e‘  (147) 

An  +  iBn  '  T°A|s=>n  B|z=li(bT 

The  superscript  notation  in  Eqs.  346  and  147  is  for  the  purpose  of 
calling  out  the  definition  of  z  being  used  in  the  evaluation. 

We  next  treat  the  indeterminate  part  of  Eq .  138  (the  A  f  0  case). 
Applying  L'llopital ' s  rule  to  the  last  term  on  the  right-hand  side  of 

Eq.  138  results  in 


4 
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„l/NpT  zM  -*{(sin  bAT)z‘'J  +  sin  [ b ( 1  -  A)t) 
A  B  sin  a)n(T/N) 


x 


2 f 7.  -  cos  a)n(T/N)] 
2Nzn-1(zn  -  cos  hT] 


z  =  Uu>n(T/N) 


(148) 


Evaluating  Eq .  148  at  z  =  l2j!wn(T/N)  produces 


G17NGT 


z-l^d/N) 


(149) 


x  |sln  bAT(cos  (unT  +  j  sin  n>nT)  +  sln  b(:-A)I|!cos  tUnC'r/N)  a-  J  s  In  wri(T/N)  -  cos  ai„(T/N)J 
(sin  wn(T/N) ]N[cos  (unT  +  j  sin  wnT  -  cos  bT] 

Simplifying  Eq .  149  gives 


.  .  ,T/N  Tl 

^TJBn  *  GA  s 


'-MUfjCT/K') 


(sin  bAT)  (ms  o.-T)  +  sin  bC  -  4)T  !  jUi...  i%T)  (sin  bur) 
N  sin  bT  ‘  . 


gT/ncI! 


Iz-lZa^WN) 


sin  bT  cos  bAT  +  j  sin  a^T  sin  bAT 
N  sin  bT 


(.T/NrT  I 


IZ-I^CT/N) 


•  cos  LAT  +  1  sin  bAT 
N 


T/N  Tj 

CA  S 


rT/NrT 
'■'A  CB 


AbT 


.L4*u„<T/N) 


Is-UiUiCT/K) 


(150) 


S“vlb 


Thus,  the  only  new  element  added  is  the  factor  evaluated  at  the 

input  frequency.  Since  A  is  generic,  we  draw  the  same  conclusion  for 
every  other  element  of  (W^R)^  and  hence  combining  Eqs.  146  and  150 
gives 
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An  +  J®n 


(gw2  )t/n 


£„sT/N 


2ce 


z=ljj:wn(T/N) 


(^2*GlW1jT 


A  cT 
Z=50O  1 

z=1j£u>tiT 

zsi<bT 


Wl*|s= 


jb 


(151) 


wn 


b 


n  = 


(152) 


As  in  Eq.  146,  letting  N  +  “  gives  the  "continuous"  result 


An  "1  Jbn 


GW  2 
“ 


S-  •’(On 


(w2*Glwl)T 

or 


z=es^  I 

Wl* 

z=l4bt  *|s=jb 

Z  =  U(OnT 


(155) 


In  Eqs.  151  and  153  it  is  important  to  observe  that  an  option  exists 
with  regard  to  (W2^G1Wj)T;  it  may  be  evaluated  at  either  l^bT  or 
l*.conT.  However,  there  is  uu  option  with  regard  to  since  it  must  be 

evaluated  at  s  =  jb.  Clearly,  the  most  efficient  procedure  is  to 
evaluate  ( ^2 i W i )  at  z  =  l.£bT  and  at  s  -  jb,  save  the  result  in  a 

polar  format,  and  then  evaluate  (GW2/N)T/lN  (or  GW0/T)  over  the  appro¬ 
priate  range  of  n. 

It  can  be  appreciated  that  the  results  of  Section  III  are  contained 
in  Eqs.  151  and  153,  since  Section  III  considers  only  the  special  case 
of  W  |  =  1,0  (a  scalar).  Also,  note  that  Wi  and  W  2  are  symbolic 

representations  that  can  represent  either  multi-rate  or  multiple-order 
sampling  formats  (e.g.,  it  is  permissible  for  W2  =  W,  etc.).  However, 
in  arriving  at  the  spectrum  of  the  continuous  variable,  uniformly  spaced 
time  samples,  in  a  T/N  time  frame,  were  used. 


63 


"  A  -.t  ,  I'ryjintf  W  ifliMi'-r  '-i-r  «  .  • 


C,  A  CLOSBD-LOOP  APPLICATION 


The  ability  to  cast  any  multi-rate  or  multiple-order  sampling  se¬ 
quence  into  a  switch  aecomposition  format  makes  the  analysis  of  closed- 

loop  systems  a  straightforward  task.  Consider  the  system  of  Fig.  23. 

T 

Proceed  as  in  the  single-rate  case.  Since  C  *  G^W^E  ,  one  must  first 
T 

solve  for  E  : 

Et  =  (WUR)T  -  (W1*G2W2)t(W2*G1W1)t  BT  (154) 


or 

ET  -  ll  +  (wuG2W2]T(W2^G1W1)Tr1(WuR)T  0  55) 


Figure  23.  Closed-Loop  System  with  Switch  Decomposition 
Therefore 


CT/N  =  (C1W1)T''N[l  +  (WuG2W2JT(W2AGiW1)T]  1  (W^R)1'  (156) 


The  coefficients  for  the  steady-state  waveform  are.  then 


An  +  JBn 


GlWi 


s=Jun 


[I  +  (wuG2U2)'r(w2*G1W1)TJ 


-1 


gpsT 


z=e 


z=UbT 
or  z=l^wnT 


s  =  jb 
(157) 
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Thu  following  numerical  example  will  help  clarify  the  specific 
details  of  Eq.  157.  (An  alternate  solution  based  or,  the  "two-rate" 
theory  of  Ref.  L  is  contained  in  Appendix  A.)  Consider  the  system  of 
Fig.  24. 

With  reference  to  Eq .  157,  Wj  =  W  =  [1,  e  anc]  =  1.0.  The 

T/N  output  equation  for  Fig.  24  is 

CT/N  =  (G1M2W)t/n  Kt  (158) 


where 


Et  =  [I  +  (WJtG2M)T(G1M2W)T)  1(WAR)T 


(159) 


Then 


CT/N  =  (G1M2W)t/N(I  +  (W^G2M)t(C1M2W)t]  1(W)tR)T  (160) 


Tlie  steady-state  coefficients  for  Eq .  160  are  given  by 


+  J'Vi 


GiM2W 


s-^JWn 


[I  +  (W*C2M)T(G1M2W)T] 


-1 


A 

z»eb  1 


z=UbT 

z--l4<j)nT 


■jb 


(161) 


Figure  24.  Closed-hoop  Numerical  Example 
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where 


=  b  + 


2Ttn 


n  =  0, 1,2, . . . ,N-1 


(162) 


The  inverse  expression  in  Eq .  161  is  formed  using  the  following 

steps  (z  =  es^): 


(W*G2M)t 


where 


! 

i 

1  -a-»T|T 

(1  -  e~T/2)z  q.  (e~T/2  —  e~T) 

t 

esT/2 

s  ( s  +  1  )  j 

z  -  e-^ 

(163) 


1  _  e  s  T 

T  1 

Z  “  1 

l 

1 

s  (s  +  1 ) 

z 

s 

s  +  1 

(164) 


and 


es'I'/2(l  -  e’"s^') 
s(s  +  1) 


T 


z  -  i  0sT/2  esT/2lT 

z  s  s  +  1  I 


(165) 


with 


>Ts 

T 

Z 

eATs 

T  —AT 

ze  n  1 

s 

z  -  1 

s  +  a 

z  -  e“aT 

(166) 


Applying  the  advanced  z-transf orm  relationships  in  Eq.  166  to  Eqs.  164 
and  165  gi  es  Eq .  163.  Next, 


(G1M2W)t 


1  -  e-sT/2 
|  s(s  +  1) 


e-sT/2j 


e-T/2  -  e-T  1  1  -  e"T/2 
_ »  _ 

z  -  e-T  i  z  ~  e~^ 

(167) 
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where 


1  _  e-sT/2 

T 

(l  -  e-sT/2)  (-1  +  e'sT/2)| 

s  ( s  +  1  ) 

s  s  +  1 

(168) 


and 


c-sT/?(i  -  e-sT/2) 

T 

e-sT/2(1  _  e-8T/2)  .  e-sT/2f +  e-sT/2)] 

s(s  +  1) 

'  1  '  T - -  -  .  ■  ■ 

s  s  +  1 

(169) 


Placing  Eqs.  168  and  169  in  an  advanced  z-transform  format  produces 


(1  -  e-sT/2)  (-1  +  e-sT/2) 

T 

-  c“l 

f_esT/2  _  esT/2 

s  s  +  1 

—  2 

t 

s  s  +  1 

s  s  +  1 


(170) 


e-BT/2(1  _  e-sT/2]  _  e-sT/2(_j  +  e~sT/2)' 

T 

_  ,  -1 

■  esT/2 

s  s  +  1 

—  z 

s 

_1 

s 


csT/2  i. 
s  +  1  s  +  1 


(171) 


Equation  166  can  now  be  applied  directly  to  Eqs.  170  and  171,  resulting 
in  Eq.  167.  Now,  assume  T  =  -ln(.81),  which  gives  a  set  of  convenient 
numbers  when  Eqs.  161  and  167  are  evaluated.  The  result  is 
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M2W)T 


[.09  j  .1] 
z  -  .81 


.0171 

!  .019 

_ 

,.009z  +  .0081 

i  . Olz  +  .009 

(z  -  .81)2 


and 


I  + 


where 


(w^G2M)T(giM2W)T  = 


z2-1.62z+.6732 
.  .009z+.0081 


.019 
z  2-1 . 6 lz+. 


■  2  -  1 . 62z  +  .6561 


6651 


I(z  -  .81)2  = 


(z  -  81)2 


0  i  (z  -  81)2 


Taking  the  Inverse  of  Eq .  173  produces 


[  I  +  (wag2m)T(g1m2w1T]-1 


z2-l .61z+.6651 


-.019 


—  ' ..009z+.0Q81  )  |  z  2_  ^  .  62y.+.6  717 


z2  -  1.61?,  +  .6822 


t-!l(z)  j  N2(z) 

n3^z)  I  N4(z) 
D(z) 
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(172) 


(173) 


(174) 


(175) 


1 T^' 


•  iititom 


Substituting  Eq  .  175  into  Eq .  161  gives 


„SV 


1  -  e-sT/2 

An  +  JBn  Ts(s  +  i ) 


c-sT/2] 


Nj(z)  ]  N?(z) 

A  oT 

z=e^  L 

Njfz)  |  N^(z) 

1 

s  =  >>n 

D(z) 

z=14bT 

ejbT/2_ 

(176) 


Performing  the  matrix  operations  reduces  Eq .  176  to  Eq.  177. 


An+ J  ®n 


where 


l-e-sT/2 
Ts  (s+1 ) 


Nj^  (z  )+N2(  z)e  z  )e  sT  /  sT/2ejbT/2 

—— 

and 


s'j'^n 

z  =  UbT 


077) 


Hi 


n 


n  =  0,1,2, ...,N-1 


(178) 


In  Eq .  177,  N^(z),  ^(z),  N^z),  N^(z),  N^(z),  and  D(z)  are  evaluated  at 
14wnT  =  14-bT  and  the  remaining  terms  at  s  =  jo>n . 

Ill  this  two-rate,  closed-loop  example  where  the  Input  sampling 
interval  Is  one-half  (T/2)  the  output  sampling  interval  (T),  two 
separate  sets  of  Bode  plots  are  required  to  present  the  magnitude  factor 
and  phase  shift  information  for  the  entire  frequency  spectrum.  That  is, 
the  fundamental  and  alias  components  are  still  directly  read  from 
individual  magnitude  and  phase  plots,  but  now  two  separate  magnitude 
plots  and  two  separate  phase  plots  exist.  One  set  of  Bode  plots  (i.e,, 
magnitude  plot  and  its  accompanying  phase  plot)  is  valid  for  a  specific 
range  of  input  excitation  frequencies  (r  =  1  sin  bt)  and  the  second  set 
for  an  alternate  range.  In  this  example,  these  two  frequency  ranges  are 
given  by 


K  tog  <  b  <  (K  +  1)ms 


(179) 
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where  U)B  =  2n/T  and  one  set  of  Bode  plots  is  applicable  for  even  K's  (0, 
2,  4,  6,  etc.)  and  the  other  for  odd  K's  (1,  3,  5,  7,  etc.).  Therefore, 
the  first  set  of  Bode  plots  is  valid  for 


0  <  b  <  uc 


2tos  <  b  <  3'jlg 


4uis  <  b  <  5ws 


(180) 


and  the  second  set  is  used  for  input  excitation  frequencies  which  fall 
in  the  following  ranges 


oig  <  b  <  2'jjs 


3ws  <  b  <  4(i)s 


5o)s  <  b  <  6ws 


(181) 


It  is  seen  that  the  ranges  covered  by  the  two  sets  of  Bode  plots 
alternate  along  the  frequency  axis. 

A  composite  magnitude  plot  for  F.q .  177  is  shown  in  Fig.  25.  The 
solid  curve  represents  one  magnitude  plot  and  the  dotted  curve  the 
second.  As  pointed  out,  these  curves  are  valid  for  the  input  excitation 
frequencies  defined  by  Eqs.  180  and  381.  Sample  printouts  are  shown  in 
Table  4  for  Eq.  177  evaluated  at  b  =  1,  40,  and  80  rad/ sec.  The  format 
of  the  printout  is  magnitude  (dB),  phase  angle  (deg),  wn  =  b  +  2x/T 
(rad),  and  an  index  number  to  indicate  the  aliased  component.  Figure  25 
was  generated  by  evaluating  Eq .  177  for  a  large  number  of  input  frequen¬ 
cies,  leading  to  the  appearance  of  a  multi-valued  Bode  plot.  However,  a 
careful  check  of  the  spectral  components  listed  in  Table  4  will  show 


TABLE  4.  DATA  FOR  COMPOSITE  MAGNITUDE  PLOT  (FIGURE  25) 


Magnitude  (dB) 

Phase  (deg) 

Frequency  (rad/sec) 
Spectral  Component  Index 


b  =  I 


-3.300 

-19.924 

1.000 

0.000 

-88.518 

■117.241 

30.818 

1.000 

-71.690 

-63.879 

-60.635 

2.000 

■107.219 

•118.466 

90.453 

3.000 

-83.587 

-64.348 

120.270 

4.000 

-116.016 

■119.718 

150.088 

5.000 

-90.582 

-64.506 

179.905 

6.000 


b  =  40 


-39.849 

-209.290 

40.000 

0.000 

-99.387 

60.648 

69.818 

1.000 

-55.701 

-210.147 

99.635 

2.000 

-116.112 

60.270 

129.453 

3.000 

63.850 

-210.362 

159.270 

4.000 

-116.694 

60.130 

189.088 

5.000 

-69.374 

-210.460 

218.905 

6.000 


b  =  80 


-51.897 

-150.751 

00.000 

0.000 

-107.245 

-428.657 

109.818 

1.000 

-61.372 

-151.057 

139.635 

2.000 

-114.780 

-428.840 

169.453 

3.000 

-67.550 

-151.180 

199.270 

4.000 

-120. C18 
-428.928 
229.088 
5.000 

-72.098 

-151.246 

258.905 

6.000 


that  the  various  aliased  components  will  alternate  between  the  top  and 
bottom  curves,  or  they  will  remain  totally  on  the  solid  curve  or  totally 
on  ttie  dotted  curve.  To  illustrate  this,  the  fundamental  and  aliased 
components  are  highlighted  with  a  star  symbol  in  Fig.  25  for  an  inpuL 
frequency  of  b  =  1.0.  The  two  magnitude  plots  in  Fig.  25  are  shown 
separately  in  Figs.  26  and  27  to  further  emphasize  that  two  completely 
independent  magnitude  plots  exist,  each  valid  for  different  ranges  of 
input  frequencies. 

This  multiplicity  in  Bode  plots  is  directly  associated  with  the 

(W  R)  term  in  Eq.  160  and  exists  for  any  general  closed-loop  system 
*  - 
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(gp)epn:jjOiui7  'SP)3Pni1la' 


with  high-i  nput:  /slow-out  put  sampling.  For  other  general  mul t i -raLo/ 
multiple  order  systems,  a  single  magnitude  and  phase  plot  covers  the 
entire  frequency  spectrum.  It  should  be  emphasized  that  the  (W  R)  term 
alone  does  not  guarantee  multiple  Bode  plots.  In  many  open-loop  and 
closed-loop  systems  the  (W  R)  ^  term  appears  and  only  a  single  set  of 

•k 

Bode  plots  is  required  ( c • g • ,  dnsc  TV  example  in  Subsection  D).  Thus, 

the  multiplicity  in  Bode  plots  is  directly  tied  to  closed-loop  systems 

with  high-input/slow-output  sampling.  For  these  systems,  the  properties 

of  the  multiple  Bode  plots  are  then  determined  by  the  (W  R)  term. 

* 

In  Eq .  161  (or  Kq .  177),  those  terms  that  are  evaluated  at  z  =  l4>'-0  T 
and  s  =  jwn  are  inherently  generic  with  respect  to  the  input  frequency. 
For  example,  given  a  sampling  freouency  of  :•)  =  2tt/T  =  10.fi  rad/sec  and 
an  input  frequency  of  b  =  1.0  rad/sec,  these  terms  contribute  to  the 
following  fundamental  and  aliased  components: 


Fundamental : 
1st  Alias: 
2nd  Alias: 
3rd  Alias: 


1.0  rad /sec 
11.0  rad/sec 
21.0  rad/sec 
31.0  rad/sec 


nth  Alias:  b  +  2nn/T  rad/sec 


Now,  for  b  =  11.0  rad/sec,  contributions  to  the  following  frequency 

components  are  obtained: 


Fundamental 
1st  Alias: 
2nd  Alias: 
3rd  Alias: 


li.fi  rad/sec 
21.fi  rad /sec 
31.fi  rad/sec 
41.0  rad/sec 


1.  +  ?nn/T  rad/sec 


74 


nth  A] 1  as : 


It  is  seen  tii.it  the  contributions  to  the  1st,  2nd,  3rd,  etc.,  aliased 
component  tor  h  =  I .0  rnd/see  identically  match  the  contributions  to  the 
tundainont.il  and  1st,  2nd,  3rd,  etc.,  aliased  component  for  b  =  li.O.  We 
could  t lion  replace  cjJn  =  b  +  2kc/T  with  simply  mn  =  b  to  obtain  the  fun¬ 
damental  and  aliased  components  for  the  Rode  plot  if  it  were  not  for  the 

T 

(W  R)  term  in  Kq  .  160-  indeed,  a  single  set  of  generic  llode  plots  can 
k 

bo  constructed  for  any  general  multi-rate/mul tiplc-order  closed-loop 
system  except  those  employing  fast-input/slow—output  sampling. 

’Hie  (W  K)  teriii,  unlike  tlie  remaining  terms  in  F.q .  160,  must  be 
evaluated  at  s  =  jb  (and  not  s  =  jw  ) .  Therefore,  this  frequency  re¬ 
sponse  term  is  solely  a  function  of  the  input  frequency  b  aid  not  a 
function  of  the  input  plus  aliased  component.  The  multiple  Bode  plots 
are  a  consequence  of'  this  singular  dependence  on  the  input  frequency  b 
and  l iui  result  ing  cyclic  or  periodic  nature  of  the  elkf/2  element  in 

! 

I  s  =  jb  t.sV2 


s-in 


(182) 


For  this  example,  this  periodicity  is  at  the  input  sampling  frequency 
[  2n  /(T/2  )  =  2i>s  ]  since 


,,  jbT/2 


•2  a/ (T/ 2) 


cos  2t  +  j  sin  2n 


(  183) 


and  the  contribution  from  the  (W  K)'  term  then  repeats  itseli  lor  input 

* 

frequencies  equal  to 


b  -  2nn/ (T/2 )  1 1  -  1,2,3,...  ( 18b) 

Tlifs  cyclic  behavior  shows  up  directly  in  the  composite  magnitude  plot 
in  Fig.  2  'j  where  the  t  undarienta  I  frequency  rang'1  covered  by  the  combi¬ 
nation  of  both  Rude  plots  equals  the1  input  samp  ling  frequency,  2‘i>f  ■» 

bV.tilb)  rad /sec .  However,  each  lk.de  plot  contains  aliased  componeuts 


7  ’> 


—  — * - - - -  *  ‘  '  Jr'*  1-  f,.-,- 


which  are  a  function  of  the  output  and  not  the  input  sampling  frequency 

(l.e.,  b  +  2n/T,  b  +  4ir/T,  b  +  6tt/T,  etc.)* 


».  GKSEBAL  RiiSUI-rS 

Some  general  c„.  usions  can  be  formed  by  extending  the  two-rate 
vector  switch  decomposition  example  from  the  previous  subsection.  The 
final  expression  (Eq.  161)  for  the  spectral  coefficients  is  repeated 
below  for  easy  reference. 


C  i  M  yW 


“  t 


s-J“a 


[I  +  (W„G2M)T(G1M2WjTj-1 


A  qT 


E=  l<wnT  */ss-jb 
a«UbT 


(185) 


Vector  sv/itch  decomposition  modeling  for  a  different  input  sampling 

interval  only  requires  alterations  to  the  W  and  W  vectors.  This  change 

* 

does  not  disturb  thoi'e  terms  evaluated  at  s  -  jwp  and  z  =  ljfmnT.  It  is 

known  that  the  W  term  evaluated  at  s  *  1b  dictates  the  structure,  of  the 
* 

frequency  response  spectrum  and  the  Bode  diagrams.  As  was  shown,  for  a a 

Input  sampling  frequency  that  is  twice  the  output  sampling  frequency, 

two  separate  sets  of  Bode  plots  are  required  as  a  result  of  the  e^T/'- 

elcment  in  the  term.  If  this  Is  extended  to  an  input  sampling 

interval  of  T/3,  the  W  vector  will  contain  the  following  elements: 

* 


w 


* 


1 

esT/3 

(.ST2/3 


086) 


Both  e-^/^  and  elements  now  appear  in  the  spectral  coefficient 

expression.  The  factor  is  periodic  in  b  =  2irn/(T/3)  and  eJki2/3 

In  b  "  2irn/(2T/3).  This  periodicity  results  in  not  two  but  three  sepa¬ 
rate  sets  ol  Bode  plots  to  cover  the  entire  frequency  response  spec¬ 
trum.  Again,  as  in  the  T/2  case,  all  alias  components  will  be  directly 
related  to  multiples  of  the  output  sampling  frequency,  <u  «»  2«/T  (l.e., 


b  +  2ir/T,  h  +  4it/T,  b  +  6ti/T,  etc.).  The  individual  range  of  input 
frequencies  that  are  applicable  to  these  three  sets  of  Bode,  plots  is 
given  by: 


(k)ii>s  <  b  < 
(k  +  l)ws  <.  b  < 
(k  -r  3/ 2)<os  <  b  < 


(k  +  l)ua 
(k  +  3/2)u>s 
(k  +  3  )(0g 


u)s  -  2r.  /  T 


k  =  0, 3, 6, 9, etc. 


(187) 


The  fundamental  frequency  range  covered  by  all  three  Bode  plots  again 
occurs  at  the  input  sampling  frequency,  3aos  • 

Further  extensions  of  this  example  provide  the  following  general 
conclusions  for  fast-input/slow-output  sampled  closed-loop  systems: 

«  The  number  of  separate  sets  of  Bode  plots  re¬ 

quired  to  display  the  entire  frequency  spectrum 
is  equal  tc  the  number  of  elements  in  the  W  | 
switch  decomposition  modeling  vector.  *  s~jb 

*»  The  frequency  range  that  is  applicable  to  each 

Bode  plot  is  determined  by  the  periodicity  of  the 
Wjs  =  jb  elements. 

<»  The  alias  components  in  all  Bode  plots  are  a 

direct  function  of  the  o  itput  .sampling  frequency, 
o>s  =  2ir/T  (i.e.,  b  +  wg ,  b  +  2<a»s  ,  b  +  3m  ,  etc.). 

*  The  fundamental  frequency  range  covered  by  the 

combination  of  Bode  plots  is  equal  to  the  input 
sampling  frequency. 

In  addition  to  theje  general  conclusions ,  the  fast-input/slow-output 
sampling  scheme  gives  rise  to  new  frequency  response  identities.  These 
new  identities,  along  with  those  presented  in  Ret.  1,  allow  the  frequen¬ 
cy  response  expressior  (A,  +  jBn)  to  bo  written  directly  from  inspection 
of  cacti  re tm  in  the  output  equation.  Any  general  multi-rate/ 

multiple-order  sampling  architecture  can  be  handled  with  this  complete 
set  of  identities.  This  includes  vector  switch  decomposition  modeling 
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of  mull: i-rate/multiple-order  systems,  computational  delays,  and  nonsyn- 
chronous  sampling. 

For  tlie  two-rate  example  in  the  previous  subsection  (Fig.  24),  the 
C^N  output  equation  using  multi-rate  theory  (Eq.  A-6),  and  the  CT//N 
equation  using  vector  switch  decomposition  modeling  (F,q  .  160)  follows: 

CT/N  *  (ciM2)T/N  rT/2  -  (GiMyJT/N  (g2m)T/2 

x  |l  +  ((c1M2)(g2m)t/2]T|  1(g1M2Rt/2)7  (188) 

CT/N  =,  (g1M2w)t-/n[  I  +  (w^C2M)T(G1tl2w)Tl_1  (w*r)t  (189) 

In  the  CT/N  equation  arrived  at  using  multi-rats  theory  ( F,q .  188),  the 
first  term  on  the  right-hand  side  introduces  spectral  components  at 
to  =  b  +  2xn/(T/2),  n .  *=  1 ,  2,  ••*,  and  the  second  term  introduces  spec¬ 
tral  components  at  =  b  -F  2nk/T,  k  =  0,  1,  2 .  The  unique  factor 

in  this  expression  is  (GjM2R^^2) ‘  which  corresponds  to  the  (W  R)2"  factor 
in  the  CT/N  equation  based  on  vector  switch  decomposition  modeling  (Eq . 
189).  These  two  factors  directly  determine  the  multiple  Rode  plot 
structure.  The  new  frequency  response  identities  for  these  two  factors 
are 

•  £,>sT/M 

(CiM2Rt/2)t  [ g'rt/M]T  G'(z)jZ  "  (190) 

|z*UbT/M 

(war)t  wj  (191) 

I  s= jb 

All  the  remaining  terms  and  factors  in  the  ^  equation  can  be 
handled  by  the  foliowing,  identities  from  Ref.  1. 
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Finite  N 


<GM)t/N 


l  (GM)t/N(z) 


zScsT/N 
z  =  Uun(T/N) 


C, 


T 

A 


A  o  T 

7.=CS  1 

z=14bT 


GT/Mk 

K 


z=esT/Mk 

z=l4o>n(T/Mk) 


n  +  -iBn  = 

TT  (CM)  T/f,(z) 

z=esT/N 

cj(z) 

h 

z-l4«n(T/N) 

A 

A  cT 

z=es 
ii  =  i4b 


x 


k 

n 


cj/rik(z) 


zftesT/Mk 

z=14ai,((T/Mk) 


K’ 


«5 


(gm)t/k  —*■  ^  (on)  (s )  I 

■s“J“n 


‘5 


CA(z) 


z=e«T 

z=UbT 


T/ 


n. 

K 


T/M, 

Gk"-(z) 


z=l4iDn(T/Mk) 


xn  + 


lBn 


i  1 

A  o  T 

•/  sn  r>  1 

Y  ( c  m  )  ( s )  I  ; 

ls=jw„J 

CA(X> 

e.  t 

z=14b T 

X 


k 


jz=osT/nk 

|z=l4«n(T/Mk) 
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n  =  nQ ,  n0  +  1,  ...  0,  1,  2,  N  -  n0  -  1 

no  =  ~^b/u)s)INT 

Ihe  notation  used  in  these  identities  indicates  the  definition  of 
"T"  to  he  used  in  calculating  each  z-transform  expression  and  the  subse¬ 
quent  evaluation  performed  to  calculate  coefficients  for  the  fundamental 
and  alias  terms  in  the  output  waveform  (An  +  jB  )«  That  is: 

•  Take  the  T/N  z-transform  of  (GM)  and  evaluate  at 
z  =  Ua>nCT/N). 

•  Take  the  T  z-transform  of  G^  and  evaluate  at 
z  =  14-bT. 

•  Take  the  T/M^  z-transform  of  each  term  and 
evaluate  at  z  =  l4^n(T/M^). 

For  the  (GM)  term,  calculate  the  normal  z-transform  expression  and 
replace  each  "T"  in  this  expression  with  "T/N."  Follow  the  same  proce¬ 
dure  for  G^  and  G^  with  "T"  replaced  with  "T/M^"  in  the  resulting 
z-transform  expression  for  the  G^  terms.  With  (GM),  GA,  and  in 
z-transform  form,  each  term  Is  then  individually  evaluated  at  the  appro¬ 
priate  value  of  "z"  as  indicated  by  the  notation. 


Example 


CT/IOO  .  (CM) T/ 100  GT  GT/10  gT/20  gT/50  rT 

A  J.  I  J 

izSesT/  J.00 

An+.Pn  (GM)T/ 100(B)  Gj(z) 


(200) 


_  A  _  «  T 

~  ti 


gT/J0(z) 


I  z=l4wn (T/100)  *'  j  z=l4bT 


A  r.T  /  1  n 

iz=eJi' 


I  z=14<*)n(T/ 10) 


IZ“esT/20 


IZ-MT/5C' 


-  Gt/20(z)  gT/50(z) 

z=l*a>n(T/20)  J  z  =  Hwn(T/50) 


(201) 
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b  +  (2Trn/T) 


(0n 

n  ~  Qq  ,  n,-,  1 ,  •  •  •  i  0,  1  ,  2 ,  •  •  .  ,  N“Dq— I 
no  ~  ~(^b/ws)^^T 

It  is  apparent  that  once  the  C^1  ^  expression  has  been  determined  for 
the  system,  simple  substitution  produces  the  exact  frequency  response 
expression  for  either  the  finite  N  case  or  the  limiting  case  of  N  ■+  «>, 
Using  this  expression,  the  coefficients  (AR  +  jBn)  for  the  fundamental 
and  its  aliases  present  in  the  output  waveform  can  then  be  determined. 
A  summary  of  the  steps  required  follows: 

•  Obtain  the  output  equation  using  block 

diagram/signal  flow  algebra. 

•  Apply  the  ident ities  outlined  in  this  subsection 
to  each  term  in  CT'N  to  obtain  the  frequency 
response  expression  +  jBn> 

•  Calculate  the  individual  z-Ltciusforms  using  tne 
appropriate  definition  of  z. 

«  Evaluate  each  term  in  the  A^  +  jBn  expression  at 
appropriate  values  of  z  and  s, 

E.  APPLICATION  TO  SIMULATION  EE808  ANALYSIS 

There  is  a  natural  extension  of  the  previous  results  which  makes  the 
switch  decompose tion/f requency  response  technique  ideal  for  the  error 
analysis  of  simulations  (either  simulations  of  continuous  systems  or  the 
simulation  of  continuous  plants  under  the  influence  of  direct  digital 
control).  One  may  obtain  a  clearer  understanding  of  this  statement  by 
considering  the  four  diagrams  of  Fig.  28. 

The  situation  described  by  Case  I  in  Fig.  28  assumes  a  controlled 

9 

element  (l/s^)  under  the  influence  of  a  continuous  feedback  controller 
with  an  (idealized)  compensation  network  2(s  +  1 )  in  the  forward  path. 
Case  li  depicts  the  same  controlled  element  under  the  influence  of  a 
discrete  feedback  controller  which  smooths  the  output  of  the  digital 


81 


computer  with  a  ZOU  (passing  on  a  "staircased"  signal  to  the  control 
point).  The  discretized  version  of  2(s  +  1)  was  computed,  using  the 
first  back  difference  algorithm  resulting  in  the  model  (42z  -  40)/z  (at 
a  sampling  rate  of  20  Hz). 

Case  Ill  depicts  a  simulation  of  Case  II  (or  I)  wherein  the  plant, 
l/s2,  is  modeled  using  a  "substitution-for-s"  rule  (Tustin  transform). 
The  output  variable  0  is  modeled  as  the  output  of  a  storage  register 
(ZOH) . 

Finally,  Oise  IV  depicts  a  common  situation  wherein  one  part  ol  a 
simulation  is  coded  for  one  computer  while  another  part  is  coded  on  a 
second  computer.  Typically,  the  computers  are  working  in  slightly  dif¬ 
ferent  frame times  and  therefore  will,  on  occasion,  pass  "old  data"  back 
and  forth.  In  this  example  it  is  assumed  that  the  compensation  is 
modeled  in  a  0.05  sec  time  frame,  while  the  plant  is  modeled  in  a  0.075 
sec  time  frame.  Data  transfer  between  the  two  computers  is  via  appro¬ 
priate  buffer  registers,  modeled  as  a  ZOH  in  a  T/3  time  frame  (M-j)  and  a 
ZOH  in  a  T/2  time  frame  (f^). 

This  completes  the  problem  description;  we  may  now  write  the  appro¬ 
priate  equations  for  each  case  and  discuss  the  analytical  difficulties. 

Case  1  is  straightforward  since 

C( s )  =  [ I  +  ClG2]“1  C1C2R  (202) 

Case  II  appears  straightforward  but  does  present  a  contradiction  in  the 
output  equation 


I  + 


(CjH)T 


QW 


(203) 


TIiat  is, 


V 


i  -  r 

s 


(z  +  1) 

2  (z  -  l)2 


(204) 
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is  readily  computed  using  a  transform  table.  However,  it  was  not  the 
intent  of  the  designer,  who  used  a  substitution  technique  to  model 
2(s  +  1)  as  (42z  -  40)/z,  to  compute 

GT  =  [2(s  +  1) ]T  (205) 

T 

using  the  z-transf orm.  Indeed,  what  does  [s]  mean?  Clearly,  the 

intent  was  to  assign  Gj  =  (42z  -  40)/z,  which  is  legitimate  if  one 

Interprets  z  *  as  a  delay  operator. 

This  difficulty  surfaces  again  in  Case  III  since 

CT  *  [l  +  GTGT]"1  GjGTRT  (206) 

after  taking  due  note  that  MT  =  1. 

Now  it  is  necessary  to  interpret  both  G?  and  as  given  functions 
of  z  rather  than  z-transtorm  operations.  For  example,  assign  (via  the 
Tustin  transform) 


GT  .  T2  (z_+  1)2 
1  *  (z  -  1)2 


(207) 


rather  than  (via  the  z-transform) 


1_IT  _  Tz 
s2  (z  -  1)2 

Difficulties  in  the  "assignment"  procedure  also 
write  the  output  equation  for  Case  IV  (see  Fig.  29). 


(208) 


surface  when  we 


Figure  29.  Case  IV,  Switch  Decomposition  Model 


C  =  (M2W2)[I  +  (W2*GlW2)T(w2*H3W3)T(W3jkG2W3)T(W3#(M2W2)T]  1 

x  (w2*GlW2)T(W2*M3W3)T(w3*G2W3)T(W3jkR)T  (209) 

wiiere 

W2  =  [ 1 »  e-s^/2J  t  W3  =  [  1 ,  e-sT/3,  e-sT2/3]  (210) 


The  meaning  of  some  terms  in  Eq.  209  are  clear;  others  are  not.  For 
example,  a  straightforward  computation  yields  (see  Appendix  B) 


(w3*M2w2)r 


1  0  0 

(  w2am3w3)t  *= 

0  1  0 


(211) 


However,  the  meaning  of  (W3*G2W3)T  or  (W?*GjW2)T  is  not  clear.  When  the 
designer  specified  G2  =  (42z  -  40 )/z,  he  intended  (for  a  T  time  frame) 
the  recursion  relationship 


X2(n)  »  42 (Rn  -  Cn)  -  40(Rn„i  -  Cn_x) 


(212) 
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) 


In  a  like  manner,  the  Tustiin  substitution  for  l/sz  was  meant  to  yield 


X4(n)  =  2X4(11  -  1)  -  X4(n  -  2)  +  yg-  1X3(11)  +  2X3(11  -  1)  +  X3(n  -  2) 


Clearly,  the  information  available  in  Case  IV  is 


(213) 


qT / 3  _  42z  -  40 

2  z 


T/3  *  0.05 


=  esT/3  (214) 


and 


GT/2 


t£  (z  +  1)2 
16  (z-  1 )- 


T/2  =  0.075 


=  esT/2 


(215) 


r 

Using  only  the  given  computer  code,  how  does  one  compute  (W^  G2W3)  ar*d 
(W ^Gi W2 ) This  problem  is  treated  in  the  next  subsection.  In  order 
to  complete  the  analysis  of  the  present  problem,  the  reader  is  asked  to 
temporarily  accept  the  assignments 


(w3*g2w3)T 


42  0 

-40  42 

0  -40 


-40/z 

0 

42 


-  ■  sT 

-  e-  1 


(216) 


and 


^2*g1^2^  = 


t2  ^ 

yj~  (z2  +  6z  +  1) 
t2 

—  z(z  +  1) 


(TZ/4)(z  +  l)/z 

li  (;2 


16 


(  Z  2  4-  gz  +  1  ) 


(z  -  i)2 


(217) 


With  these  assignments,  all  terms  in  Eq.  209  are  defined  and  the  Bode 
plots  descriptive  of  Cases  1-1V  can  be  computed.  The  results  (magnitude 
plots  only)  are  shown  In  Figs.  30  through  35. 
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Inspection  cr  Fig.  30  discloses  no  surprises.  The  digitally  con¬ 
trolled  system  Is  a  reasonably  faithful  reproduction  of  the  analog 
system  until  the  folding  frequency  (approximately  equal  to  62.8  rad/sec) 
is  passed.  Notice  that  in  the  discretely  controlled  system,  minimum 
response  points  in  the  Bode  plot  (notches)  occur  at  multiples  of  sam¬ 
pling  frequency  (approximately  equal  to  125.66  rad/sec). 

The  comparison  of  the  single-rate  simulation  against  the  baseline 
design  exhibits  fidelity  over  a  shorter  low-frequency  range  (Fig.  31). 
Of  particular  Importance  is  that  the  aliased  bands  exhibit  a  much  higher 
amplitude  response  than  did  the  aliased  bands  of  the  digitally  con¬ 
trolled  system  (Case  II).  Moreover,  the  "notches"  now  occur  at  multi¬ 
ples  of  the  folding  frequency  rather  than  the  sampling  frequency.  In 
addition,  there  are  very  sharp  notches  which  occur  close  to  odd  multi¬ 
ples  of  the  folding  frequency;  these  are  a  consequence  of  the  zeros  of 
the  Tustin  transform  Introduced  by  (z  +  1)^.  A  direct  comparison  be¬ 
tween  Cases  II  and  III  is  given  in  Fig.  32. 

Figure  33  compares  the.  two  independent  processor  case  (Case  IV) 
against  the  continuous  baseline  design.  Large,  sharp  resonant  peaks 
have  been  introduced  in  the  aliased  bands  and,  in  addition,  there  is  a 
significant  overshoot  in  the  first  fold.  A  comparison  between  the  two 
rate  simulation  (Case  IV)  and  the  digitally  controlled  continuous  system 
(Case  II)  is  given  in  Fig.  34.  Figure  3b  (a  very  "busy"  figure)  com¬ 
pares  all  the  cases. 

There  are  significant  differences  in  the  spectral  content  of  the 
four  cases  which  would  be  hidden  if  one  only  looked  at  the  sampled  spec¬ 
trum  (that  is,  looked  only  at  the  frequency  content  from  zero  to  the 
first  folding  frequency).  Even  in  the  first  fold,  there  is  a  signifi¬ 
cant  difference  in  the  Bode  plot  of  the  continuous  case  and  the  two-rate 
simulation;  the  reason  for  the  added  overshoot  in  the  two-rate  simula¬ 
tion  will  be  discussed  in  the  next  subsection.  After  this,  we  return  to 
the  theoretical  issue  of  computing  such  terms  as  G2W3)  w*ien  only  a 

assignment  in  z  has  been  given  (or,  equivalently,  the  recursion 

2 

equation  is  given). 
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•ase 


Cowpar ison ;  Case  II  and  Case  III 


Figure  33.  Magnitude  Plot;  Case  I,  Case  IV 


92 


Figure  34.  Comparison,  Case  II  and  Case  IV 


iparison.  Ail  Cases 


F.  AN  IMPORTANT  AI.GORTTHM 


Of  fundamental  concern  is  the  evaluation  of  terms  such  as  (W  GW)^ 
when  G  is  specified  only  in  terms  of  software  statements  (i.e.,  given  as 
a  function  of  z)  .  As  a  first  step,  consider  how  one  might  evaluate  the 
system  of  Fig.  36.  As  indicated,  it  is  permissible  to  insert 


Figure  36.  A  Phantom  Sampler  Formulation  of  a 
T/N,  T/M  Sampling  Format 

a  phantom  sampler  between  G  and  the  T/M  sampler  since  T/MN  and  T/M  are 
integerly  related.  Thus, 

CT/M  =  [r,RT/NjT/M  _  f  CT/MN  rT/N)T/H  (2l8) 

M  and  N  are  integers. 

To  compute  the  entries  of  Eq.  21fi  one  may: 

1)  Compute  R  by  table  lookup. 

2)  Replace  T  by  T/N, 

3)  Replace  z  (in  R)  with  zM. 

4)  Compute  G  by  table  lookup. 

5)  Replace  T  with  T/MN. 

6)  z  remains  the  same  (in  G). 

7)  Compute  using  Eq .  219. 

cT/M  '  27T  f  G(>,)  r(PM)  — Vi  ^  (219) 

J  J p  z  -  p H  P 

Equation  219  is  mechanized  in  a  digital  computer  program  (TXCONV)  which 
is  discussed  in  Volume  II. 
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For  example,  let 


z 

G 


R 


esT/6 

1 

s  +  2 

1 

s  4  I 


M  =  3 

gt/6 

RT/2 


N  *  2 


z 


(220) 

(221) 

(222) 


so  that 

CT/3  *  J_  f _ _ _ _ & 

2*j  J  (p  -  e"T/3)(p3  -  e-T/2)(z  -  p2)  p 


1  f _ _ g_^z _ _  , 

2"J  J  (P-e-T/3)(p-e-T/6j(p2+e-T/6p+e-2T/6)(z_p2)  ap 


=  res 


j  pa=0~X/  3 


+  res  I 


lp=e-T/6 


+  res 


p—(e  T/6)/2+je~T/6(/3/2) 


+  res  I 


P=-(e“T^,/2)- je  ^//^>(»'3/2) 


(223) 


Thus,  there  are  tour  residues  to  evaluate.  Relegating  the  numerical 
details  to  Appendix  C,  we  proceed  directly  to  the  final  result, 


_ z2(z2+e-5T/6) _ 

(  z-e“-^/3)(Z-e~3’/3  j^z2+e-T/324.e-2T/3l 


z  =  esT/3  (224) 


and  observe  that  the  T/2,  T/3  sampling  format  has  produced  additional 
lightly  damped  modes  in  the  output  response.  The  reality  of  the  addi¬ 
tional  modes  can  be  better  appreciated  by  first  plotting  the  continuous 
variable  C(t)  and  then  picking  off  the  T/3  sample  points.  This  is  done 
in  Fig.  37.  Joining  the  sample  points  with  a  smooth  curve  emphasizes 
the  lightly  damped  nature  of  the  response.  This  effect  was  also  present 
in  the  two-rate  simulation  analysis  of  the  previous  subsection  (recall 
the  additional  "overshoot"  in  the  first  fold  of  Fig.  33). 
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G.  COMPUTATION  OF  (W  GW)T,  GT/N  GIVEN 

* 

With  the  aid  of  F.q .  218,  wo  are  in  a  position  to  evaluate  configura¬ 
tions  sucli  as  those  of  Fig.  38  (the  vector  switch  decomposition  equiva¬ 
lent  is  shown  in  Fig.  39). 


Figure  38.  Example  Two-Rate  Open-Loop  System 


Figure  39.  Switch  Decomposition  Formulation  for  Fig.  38 


From  Fig.  39, 

C  =  M2W2(W2^G2W2)'1'  ( ) T( '^3^0^113) T  (w3*R)^  (225) 

T 

As  noted  previously,  the  evaluation  of  terms  such  as  (W2  h^'3)  is 
routine  (Appendix  B) ;  our  concern  is  with  terms  like  (W^^GjWj)  when 
Gj/3  is  given  as  a  software  specification  (e.g.,  computer  algorithm). 
Focusing  on  (w3aGjWj)3,  we  find 


(226) 


(227) 


Observe  the  super-diagonal  terms.  They  are  simply  appropriate  sub¬ 
diagonal  terms  multiplied  by  e-s^.  Indeed,  it  is  only  necessary  to 


compute  the  first  column  entries  in  order  to  completely  determine 


(w3  giW3) 


Next,  observe  the  problem  statement  stipulates  Gj 


through  a  substitution  rule.  If  e“ 


is  defined  as  z 


perhaps  defined 
we  have,  using 


Eq.  218, 
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(w3ytn1w3)T  = 


gT/3  z-3(z2GT/3) 


zG|/3 


z2GT/3 


gJ/3 


zG[/3 


e~3(zg|’/3J 

z-3(z2G[/3) 


GT/3 

1 


=  pST/3 


(228) 


Although  each  term  in  Eq .  228  can  be  evaluated,  using  Eq .  219,  only  the 
first  column  entries  need  be  computed  (i.e.,  the  three  entries  of  column 
one  are  sufficient  to  define  all  nine  entries). 

To  summarize, 


(w3*g1w3)T  =  tw3*Gl'/3w3) 


(229) 


For  example,  recall  the  G|’'/3  transfer  function  from  the  simulation  exam¬ 
ple  (Case  IV)  had  the  form 


gT/3  =  a°Z  +  31 

'1  z 


z  =  esT/3  (230) 


One  must  therefore  compute 


/  anz  +  at 

\  z  / 


/z(a0z+aj)yr  {z2[hqz  +  ai)\T 
\  Z  /  ’  \  z  / 


(231 ) 


in  order  to  define  the  3x3  matrix  (W^GjW^) 


To  illustrate, 


^a0z  +  a1jT 


_1__  C  aQP  +  al  #  z  dp 
2,j  ir  P  z  -  p3  P 


d_  /(aQP  +  al)z\ 
dp  \  z  "  P3  \ p=0 


a0 


(232) 
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U0z  +  ai)T  .  —  f 


l  r  (aop  +  ai)z 

2lTj  Jp  p(z  -  p3j 

(a0p  +  a^zl 

- 3 -  =  al 

2  "  P  p=0 


[  (z ) (agz  +  ai)]T  =  ~  J 


1  f  UOP  +  al)2 


z  -  p  ^ 


(no  poles  enclosed).  Therefore, 


(w3*g2w3)t 


ao  0  aj/z 

al  a0  0 

0  ai  ao 


=  PsT 


Setting  ag  =  42,  =  -40  verifies  Eq.  216  of  the  previous  subsection. 

The  other  assignment  made  was 


(T2/16)(z2  +  6z  +  1)  [ (T2/4) (z  +  1 ) ] / z 


j  (T2/4)[z(z  +  1)]  (T2/16)(z2  +  Gz  +  1  )J 

(W2*G1W2)T  = - ; - — 3 - (236) 

(z  -  l)4 
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I 


*»T/2  »  JL_ 


T2  (z  +  1); 


(z  -  1); 


z  =  esT/2 


(237) 


To  verify  this,  write 
.T/2„>T  _ 


("2*G1  w) 


"  ■ 

1 

H  *->l  Cj/2 

T  (t2  (z  +  l)2 

1  Z  ^ 

! 

)  16  (z  _  1)2 

z  1 

z 

( 

J 

The  first,  terra  of  column  one  is 


16(z  -  l)2j 


The  second  term  of  the  first  column  is  easier  to  compute: 


I  j  T^_  p(p  +  l)2 


z  dp 

-  — -  =  res 


2irj  J  16  (p  -  1)2  (z  -  p2)  p 


3=1 


T2  d_  z(p  +  l)2 
16  dp  (Z  _  p2) 


T^_  z(z  +  1) 
4  (z  -  l)2 


P=3 


(240) 
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'  ***-.»*  -  4  •.>  &H.  *  ■  1  •*» 


pS 


l 

r  t2  (p+d2 

z 

= 

res 

res  1 

i  '  ■ 

2rrj 

Jr  16  (p_i)2  j 

;-p2 

P 

p— 0  1 

p=1  ! 

T2  . 

T2  d  (p+l)2z 

T2 

1  4*  - 

8z 

i 

i 

16  T 

16  dp  p(z_p2) 

p=l 

”  16 

1  T 

1 

(^-D2. 

! 

j 

T2 

(z  2  +•  6z  +  1) 

z  =  esT 

(239)  ; 

16 

(7.  -  l)2 

,  1 

I  ! 

The  only 

remaining  terms 

are 

the 

2,2  term  (which  is 

the 

same  as 

the 

1,1 

term)  and 

the  1,2  term. 

The 

1,2 

term  is  the  product 

of  z 

with 

the 

2,1 

term.  Thus,  we  verify 

Eq. 

236 

of  this  subsection 

and 

Eq.  217 

of 

the 

previous  subsection. 

H.  SECTION  SUMMARY 

Important  extensions  to  the  frequency  response  concept  for  single¬ 
rate  systems  have  been  given.  These  extensions  make  it  possible  to 
investigate  the  spectral  content  of  the  continuous  variables  of  multi¬ 
rate  and/or  multiple-order  configured  systems.  In  addition,  an  impor¬ 
tant  relationship  was  deduced  which  makes  it  feasible  to  apply  the 
switch  dec jmposition/f requency  response  technique  to  the  error  analysis 
of  simulations  of  continuous  systems  and/or  discretely  controlled  con¬ 
tinuous  systems.  An  illustrative  example  was  used  to  demonstrate  the 
significant  spectral  differences  that  occur  when  a  closed-loop  system 
(either  an  analog  or  digital  controller)  is  compared  with  an  all-digital 
simulation  of  the  closed-loop  system.  The  example  also  treated,  as  a 
particular  case,  the  problems  encountered  when  a  simulation  software  is 
put  up  on  two  different  (independent)  computers,  each  working  in  a 
slightly  different  frame  time. 

At  this  juncture,  the  important  spectral  analysis  tools  have  been 
developed,  and  we  return  (in  the  next  section)  to  the  study  of  the 
multiple-order  vector  concept. 


102 


SECTION  V 


A  PSEUDO  MEASUREMENT  EXERCISE 


The  pseudo  measurement  approach  achieves  certain  well-defined  goalB, 
as  demonstrated  In  Section  II.  This  section  addresses  the  Important 
Issue  of  whether  or  not  the  gains  required  will  be  reasonable  when  plant 
dynamics  typical  of  modern  aircraft  are  used. 

To  provide  a  digital  controller  baseline  for  the  numerical  ranges 
involved,  we  first  "emulate"  an  analog  controller  for  the  YF-16.  Next, 
the  pseudo  measurement  approach  provides  alternate  loops  for  the  basic 
airframe.  The  numerical  range  of  the  compensation  can  then  be  meaning¬ 
fully  evaluated,  though  the  Intent  is  not  to  force  the  pseudo  measure¬ 
ment  design  to  "copy"  the  emulated  baseline  controller.  A  basic  frame 
time  of  T  =>  0.04  sec  (25  Hz)  Is  used  throughout. 

A.  ANALOG  CONTROLLER - YF-16 


An  analog  flight  controller  for  the  longitudinal  YF-16  is  given  In 
Fig.  40  (M  =  0V6,  h  =  .30,000  ft).  This  figure  (and  the  transfer  func¬ 
tions  for  the  aircraft  model)  were  developed  using  data  from  Ref.  5. 

A  listing  of  transfer  functions  for  this  flight  condition  can  be 
found  in  Appendix  D. 

A  digital  controller  configuration,  using  a  ZOH  together  with  com¬ 
pensation  defined  by  the  Tostln  transform,  is  shown  in  Fig.  41. 


For  frequency  response  purposes,  note  for  example 


/  N<5  \  T  T 

«  =  (—  CiMjG26ec 


(241) 


where 


T  T 
GAR 


(242) 
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That;  is,  the  coefficients  for  the  frequency  response  components  are 
defined  (for  a): 

(G2cJ)(z>|  (243) 

S=ju)n  |z  =  l  (0nT 

T 

There  is  a  need  to  know  the  Tustin  transform  of  such  elements  as  , 
T  T 

G3 ,  G4 ,  etc.  A  significant  computational  burden  can  be  neatly  side¬ 
stepped  at  this  point  by  noting  that  the  bilinear  w'-transformation 
applied  to  a  z— domain  transfer  function  (which  itself  was  derived  using 
the  Tustin  transform)  yields  the  same  form  as  the  s-domain  function. 
For  example,  the  compensation  network 


An  +  JBn 


»6, 


“GiM 


- 1 . 2s  (s  +  5)  ,  -1.2w'(w'  +  5) 

(s  +  l)(s  +  15)  (w'  +  l)(w'  +  15) 


(244) 


-10 
s  +  10 


=> 


-10 


w'  +  10 


,  etc. 


(245) 


By  working  in  w* ,  the  computational  burden  imposed  by  poles  and  zeros 
very  close  to  the  unit  circle  is  avoided.  Of  course,  Eq.  243  must  be 
modified  to  indicate  the  dependence  on  w'  rather  than  z: 


N«e  ' 

An  +  JBn  =  ~ n  G1M 


(gM)(w') 


s  — 


(246) 


w'=j(2/T)  tan  mn(T/2) 


X 

In  addition,  the  transfer  functions  [ (N6e/D)G1M]  >  etc.,  which 
define  the  effects  of  the  ZOH  on  the  open-loop  plant  must  be  tabulated. 
These  are  also  listed  in  Appendix  D. 
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B.  FREQUENCY  RESPONSE 

Using  Eqs.  243  and  246,  the  frequency  response  of  the  continuous 
controller  can  be  contrasted  against  the  discrete  controller  results 
(see  Fig.  42)  for  the  a  component  of  the  state  vector.  The  magnitude 
plot  for  the  discrete  controller  is  a  good  replica  out  to  about  70  rad/ 
sec;  however,  the  phase  angle  plots  begin  to  deviate  much  earlier  (about 
7  rad/sec).  These  results  are  typical  of  systems  designed  via  the 
emulation  approach. 


Figure  42.  Bode  Plot,  Analog  Compared  with 
Emulated  Controller 
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Phase (deg) 


Amplitude  (db) 


C.  A  COMMENT  ON  A  SIMULATION  ARTIFACT 


Suppose  the  system  of  Fig-  41  is  simulated  in  a  particular  manner. 
Specifically,  the  compensation  is  Tustin  transformed,  the  aircraft  is 
simulated  on  an  analog  computer,  but  a  discrete  model  of  the  actuator  is 
used  (Tustin  transformed).  The  Bode  plot  for  the  a  response  (Fig.  43) 
discloses  additional  "notches"  (compare  with  Fig.  42)  at  odd  multiples 
of  the  folding  frequency  (62.83  rad/sec).  These  additional  notches  are 
Introduced  as  a  consequence  of  the  algorithm  used  to  form  the  actuator 
model.  That  is,  the  Tustin  transform  of 


20  (2/T)  +  20 

s  +  20  '  20  -  (2/T) 

z  1  20  +  (2/T) 


(247) 


Figure  43.  Simulation  of  Actuator  via  Tustin  Transform 


Phose(deg) 


introduces,  because  of  the  (z  +  L)  numerator,  nulls  at  odd  multiples  of 
n/T.  This  points  up  a  general  result:  low-frequency  checks  on  computa¬ 
tional  algorithms  (perhaps  with  the  aid  of  the  "discrete"  frequency 
response)  overlook  nulls  introduced  in  the  higher  folds. 

D.  THE  PSEUDO  MEASUREMENT  DESIGN 

Assume  the  availability  of  only  an  angle  of  attack  (a)  measurement. 
A  set  of  loops  around  the  airframe  (ignoring  the  washout  network  of  the 
baseline  configuration)  is  shown  in  Fig.  44.  The  K£  feedforward  gain  is 
used,  in  the  absence  of  the  washout,  to  set  the  steady-state  gain. 
First,  survey  the  w'-plane  transfer  functions  associated  with  the  time- 


Figure  44.  Alternate  Configuration 

T 

advanced  elements  of  (W^GGjM)  (these  are  listed  in  Appendix  D) .  This 
listing  places  in  clear  view  the  w'-plane  property  of  generating  trans¬ 
fer  functions  where  the  order  of  the  numerator  and  denominator  are 
equal.  For  example,  the  transfer  function  for  a  time  advance  of  0.5T  is 
given  in  Eq.  248. 
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The  expected  result  of  having  the  phugoid  cancel  in  the  a-related 
transfer  functions,  for  a  zero  time  advance,  is  seen  to  persist  for  non¬ 
zero  time  advances  as  well.  The  greatest  effect,  as  A  varies,  appears 
in  the  high-frequency  modes. 

Comparing,  for  example,  the  transfer  function  for  A  =  0.2  and 
A  =  0.8  exposes  what  appears  to  be  an  interesting  general  trend 
(Eq.  249).  That  is,  the  numerator  polynomials  have  essentially  the  same 
coefficients  up  to  the  w'J  term,  which  suggests  subtracting  the  two  sig¬ 
nals  [  e-2sT^N^e/DjGj^M]T  and  [e'SsT^Ng^/DjGjMjT  to  form  an  effective 
transfer  function  with  a  free  w'  (see  Fig.  45). 


G.2 


C.8 


(.8525x10-7)|  (w'-K0Q8)2+(. 07455 )2] (w'-bl .82 ) (w'+50. 31 ) (w'+6 . 3x 10? ) 

D 


.  8525* 10-7w'5+.5372x10-1w'4-.O8056s3-140.1w'2-2.242w'-. 7873 

D 


(  ? Ada l 


.  4629*  1Q-3w'5+.oi  141w'4-3 ,442w'  3-14Q.  1w'2-2.261w'-.  7873 

D  (249b) 


From  Fig.  45, 


U?  =  (l+KTK^GG^jT+KTKTz-lf^e^sTGG^jT^e^sTGG^)1]}  (250) 

The  presence  of  z-^  premultiplying  the  pseudo  differentiator  term 
requires  the  modification  of  the  and  G  g  transfer  functions  by 

e-^  “  (-w'  +  50)/(w'  +  50)  since  the  synthesis  is  being  performed  in  w' 
rather  than  the  z-domain. 
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Using  these  gains  as  a  guide,  the  exact  equation  for  A(w'), 

A(w')  =  D}I  +  rTk^gG^JT  +  g|KT  [(e-2sTGG1M)T  -  ( e -8sTgg  {  m)  TJ  ( 

(253) 

is  used  to  fine  tune  the  gains.  This  results  in 


K?kT  =  -1.5 


kTkT  =  2.2 


(254) 


with  (ignoring  the  phugoid) 


A(w')  =  1.01062033W'3  +  I7.89902448w'2  +  86.32467180w'  +  165.4746084 


=  1 . 01 [ (w'  +  3.086)2  +  (2. 159)2] (w"  +  11.537) 


(255) 


The  design  appears  to  be  reasonably  insensitive  to  gain  variations 
(see  Eqs.  256-258): 


A (w ' )  =  [ (w'  +  1.36)2  +  ( 1 . 2 )2 ] (w  '  +  16.2) 


T  T 

KjKf  =  -0.7 

,.T..T  .  . 

=  t»u 


(256) 


A (w ' )  «  [(w'  +  1.003)2  +  ( 1 .45)2  J (w'  +  17,2) 


T  T 

KlKf  =  -0./ 
T  T 

K2Kf  =  0.7 


(257) 


A(w')  =  [(w'  +  6.16)2  +  (8.74)2] (  w '  +  3. 214) 


T  T 

KiKf  =  -3.0 
T  T 

K2Kf  =  4.4 


(258) 
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The  gains  used  are  comparable  in  magnitude  to  those  in  Che  baseline 
design  (refer  to  Fig.  41).  The  (r')  characteristic  equation  for  the 
configuration  of  Fig.  41  is  Ac: 

Ac  -  [(v'  +  2.154)2  +  (3.503)2] [ (w'  +  .5)2  +  (,303)2](w'  +  .00172) 

X  (w'  +  I0.7)[(w'  +  14.12)2  +  (14.48)2] (w'  +  13.55)  (259) 


The  design  which  uses  q  and  ^  feedback  has  introduced  additional  damp¬ 
ing  in  the  phugoid  term,  which  a  feedback  alone  cannot  do.  Again,  no 
attempt  to  "copy"  the  baseline  design  was  made;  the  sole  objective  was 
to  establish  feasibility  and  determine  if  reasonable  gains  were  re¬ 
quired.  It  appears  that  the  concept  is  sound  and  warrants  further 
Investigation.  Further,  the  ability  to  generate  a  "pseudo  differentia¬ 
tor"  by  differencing  time  delay  components  of  a  signal  is  interesting  in 
itself. 

The  Bode  plot  of  Fig.  46  was  generated  on  the  assumption  of  a  sinu¬ 
soidal  input  for  (see  Fig.  45). 

K.  SECTION  SUMMARY 

The  pseudo  measurement  controller,  which  gives  closed-loop  dynamics 
similar  to  an  emulated  baseline  design,  was  shown  to  be  feasible  using 
gains  of  reasonable  magnitude.  Moreover,  by  differencing  two  time 
delayed  signals,  it  was  possible  to  generate  a  new  transfer  function 
containing  a  free  w'  in  the  numerator,  making  it  a  relatively  simple 
task  to  Independently  control  two  coefficients  in  the  closed-loop  char¬ 
acteristic  equation. 

As  an  aside,  a  simulation  artifact  introduced  in  the  higher  folds, 
as  a  consequence  of  using  the  Tuscin  transform,  was  pointed  out. 
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SECTION  VI 


SUMMARY  AND  CONCLUSIONS 


The  ability  of  the  digital  computer  to  perform  a  variety  of  uisusual 
operations  on  input  data  lends  credence  to  the  idea  of  achieving  closed- 
loop  systems  with  unusual  attributes.  This  is  the  key  point  of  the 
present  study  where  the  logic  of  the  computer  is  used  to  take  a  scalar 
continuous  output  signal,  sample  it  at  a  relatively  high  rate,  and  then 
sort  it  out  (in  a  particular  fashion)  into  a  vector  of  "pseudo"  measure¬ 
ments  . 

The  concept  is  based  on  using  information  contained  i.i  a  scalar 
measurement,  sampled  at  a  high  rate,  to  generate  a  lower  rate  measure¬ 
ment  vector.  This  vector  is  likened  to  pseudo  states  since  it  effec¬ 
tively  provides  distinctively  different  information  in  the  form  of 
separate  sampling  sequences  shifted  in  time.  Vector  switch  decomposi¬ 
tion  techniques,  which  provide  a  mathematical  model  of  the  pseudo 
measurement  vector,  were  reviewed  and  extended  in  Section  II.  The 
pseudo  measurement  concept  was  defined  and  Illustrated  with  an  example. 
The  example  demonstrated  potential  since: 

*  Unity  steady-state  gain  without  forward  loop  com¬ 
pensation  was  achieved. 

•  A  five  component  multiple-order  sample  vector  was 
generated  from  a  scalar  analog  signal. 

9  Three  poles  were  assigned  using  only  one  availa¬ 
ble  analog  measurement  and  therefore  accomplished 
pole  placement  with  multiple  samples  of  a  single 
output  variable  (In  distinction  to  single  samples 
of  several  output  variables). 

The  illustrative  example  of  Section  II  did  not  address  the  practi¬ 
cality  of  the  approach,  since  the  magnitude  of  key  parameters,  such  as 
feedback  gains  (when  plant  dynamics,  typical  of  modern  aircraft,  are 
used),  was  Ignored.  This  deficiency  was  remedied  in  Section  V,  where  a 
YF-16  case  study  example  was  used  to  compare  pseudo  measurement  gains 
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against  conventional  YF-16  analog  controller  gains.  The  pseudo  measure¬ 
ment  gains  were  found  to  be  comparable  with  the  conventional  gains, 
indicating  the  approach  does  not  require  control  authority  more  exces¬ 
sive  than  Lhe  analog  controller  design. 

In  Sections  III  and  IV  the  basic  tools  needed  to  assess  the  pseudo 
measurement  concept  were  reviewed  and/or  enhanced.  The  concept  of  the. 
frequency  response  of  the  continuous  variables  of  single-rate  digitally 
controlled  systems  was  reviewed  in  Section  III  and  then,  in  Section  .  , 
given  sufficient  generality  to  treat  multi-rate  systems  (including  the 
pseudo  measurement  concept). 

An  unanticipated  dividend  accrued  from  the  analytical  development  of 
Section  IV.  As  a  "byproduct,"  an  important  relationship  was  deduced 
which  made  it  feasible  to  apply  the  switch  decomposition/frequency 
response  technique  to  the  error  analysis  of  simulations  of  continuous 
systems  and/or  discretely  controlled  continuous  systems.  An  illustra¬ 
tive  example  demonstrated  the  significant  spectral  differences  between  a 
closed-loop  system  (employing  an  analog  or  digital  controller)  and  an 
all-digital  simulation  of  the  closed-loop  system.  The  example  also 
treated  problems  encountered  when  simulation  software  is  implemented  on 
two  different  (independent)  computers,  each  working  in  a  slightly 
different  frame  time. 

Thus  the  tools  of  Section  IV  can  be  profitably  applied  to  the 
"error"  analysis  of  simulations  -  perhaps  to  predict  expected  differ¬ 
ences  between  actual  flight  test  and  man-in-the-loop  moving-base  simula¬ 
tions  . 
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APPENDIX  A 


FREQUENCY  RESPONSE  OF  FAST-INPUT/ SLOW-OUTPUT 
CLOSED-LOOP  SYSTEM  USING  MULTI-RATE  THEORY 


The  following  solution  of  the  closed-loop  system  in  Fig.  A-l  is 
based  on  the  multi-rate  theory  presented  in  Ref.  1  This  same  system  was 
treated  in  Section  IV  using  vector  switch  decomposition  modeling.  The 
objective  is  to  calculate  the  frequency  response  expression  for  the 
output  variable  The  procedure  for  this  example  is  typical  for 

closed-loop  systems  employing  fast-input/slow-output  sampling.  - - - 


Figure  A-l.  Fast-Input/Slow-Output  Closed-Loop  System 


CLOSED-LOOP  MULTI-RATE  CONFIGURATION 

The  T/N  output  equation  for  Fig.  A-l  is  given  by 

CT/N  =  (GlM2)T/N  eT/2  (A-l) 

First,  solve  for  the  ET/^  signal  at  the  input  of  the  M2  data  hold. 

ET/2  =  rT/2  _  (g2M)T/2  (GlM2bT/2)T  (a-2) 


A-l 


Next,  premultiply  Eq.  A-2  by  GjM2  and  take  the  T  transform  of  both  sides 
of  the  resulting  equation  (or  sample  both  sides  of  the  equation  at  a  T 
interval).  The  result  is 


(g1M2Et/2)T  =  [G1M2Rt/2)T  -  [(C1M2)(G2M)t/2]t  (g1M2Et/2)T  (A-3) 

Rearranging  Eq.  A-3, 

(g1m2et/2)t  =  ji  +  [  (g  iM2)  (g2m)  T/2]^  1  [gjm2rT/,2]t  (A- 4) 

and  substituting  Eq.  A-4  into  Eq.  A-2  produces  Eq.  A-3. 

eT/2  ~  -rT/2 (g2m)t/2^I  +  [(G1M2)(g2m]t/2]T|  1  (G1M2RTy'2)1  (A-5) 

Finally,  substituting  Eq.  A-5  into  Eq.  A-l  gives  the  output  equation 
Cf/,N  as  a  function  of  the  input  RT^2. 


CT/N  *  (u1M2)t/NrT/2  -  (G1M2)T/N(G2M]T/2|I+[(GlM2)(G2MjT/2jT|-1rGiM2RT/2)T 

(A-6) 


Tem  L  in  Eq.  A-6  introduces  spectral  components  at 

<«>„  =  b  +  2irn/(T/2)  ,  n  =  0,  1,  2,  3,  . . .  (A-7) 

and  Terra  2  introduces  components  at 

0Jk  =  b  +  2irk/T  ,  k  =  0,  1,  2,  3,  ...  (A-8) 

Assume  f  =  -ln(0.81);  this  gives  a  set  of  convenient  numbers  when 
Eq.  A-6  is  evaluated.  We  first  evaluate  the  terras  in  the  inverse  ex¬ 
pression 


[(G1M2)(G2M)T/2]T 


<A-9) 


Let  z  =  esT,(,2j  and  introduce  a  phantom  T/2  sampler  to  Eq.  A-9. 


Ug1M2)[g2M)T/2]T  =  [(GlM2)T/2(G2M)T/2]T 

This  mathematical  operation  is  depicted  in  Fig.  A-2. 


(A- 10) 


This  step  is  conceptually  correct  since  the  T  output  sampler  simply 
rejects  all  the  unwanted  samples  from  the  T/2  sampler.  Therefore, 


(C1M2)r/-''  = 


L  -  e-T/2 
z  - 


and 


[[GlM2)T/2(G2M)T/2] 


(  i  -  e~T/2  j 2 (  + 

z(z“  e~2’/,2  J  2 


=  p-sT/2 


z  =  e 


( A- 1 2  ) 


Equation  A-12  can  be  solved  by  calculating  the  residues  of  the  following 
expression  (see  Section  IV): 


1  f  (l  -  e~  ^  /  2  ) 2  ( p  +  j )  z  dp 

2lTj  “r  p(p  -  eT,/2  )2(z  ~  P2j  p 


(  A- 1 3  ) 


A- 3 


The  residues  for  Che  double  poles  at  p 


0  and  p  =  e 


-T/2 


are 


lp=o 


d_ 

dp 


(p  +  1  )z 


LP 


e~T/2j^2  .  p2j 


p=0 


e-T  +  2e"T/2 


-2T 


res  I 


.  A 

(p  +  l)z 

|p=e-T/2  dp 

p2(z  -  P2) 

P“e 


-T/2 


-(e-T  +  2e~T/2)z2  +  [ 3e-2T  +  4e~3T/2]z 
e-2T(z  _  e-Tj2 


Combining  Eqs.  A-14  and  A-15  gives  the  T  transform  required 


(1  -  e-T/2): 


res  1  +  res  . 

p=0  p=e-i'2 


_I _ -  e  T/ 2]  z  +  [  e  T 

e-T)2 


Substituting  Eq.  A-16  into  the  inverse  expression  in  Eq.  A-6 
for  the  inverse  results  in 


|i  +  [(glm2)Cg2m)t/2]t}  1 


_ (z  -  e~T) 2 

z2  +  [ -2e-^+( l-e~T/2] z  +  [  2e 


(A-14) 


(A-15) 

in  Eq.  A-12. 

*  2e-T/2)l 

(A-16) 

and  solving 


-2T_3e-T+2e-T/2] 
(A-17 ) 


A- 4 


For  I  =  -ln(0.81)  =  0.210721031  sec,  Eq .  A-17  becomes 


{l  +  [(gim2)[g2m)t/2]t}  1 


(z  -  0.81)2 _ 

z2  -  1 . 61z  +  0.6822 


z  =  esT  (A-18) 


Next,  calculate  the  forcing  function 

[(g1M2]kt/2JT  =  [(g1M2)T/2RT/2]T  (A-l 9) 

Equation  A-19  can  also  be  solved  using  the  residue  method  employed  on 
Eq.  A-12.  However,  a  second  method  (Ref.  1)  Is  to  use  partial  fraction 
expansion  on  the  T/2  z-transform  product  (z  =  esT/^)  in  Eq.  A-19;  take 
the  inverse  of  each  term  back  into  the  continuous  time  domain;  and  then 
resample  back  to  the  z-domain  for  a  T  sampling  interval  (z  =  es^>.  In 
practice,  the  intermediate  step  of  computing  the  generating  time  func¬ 
tion  is  merely  a  convenient  conceptual  step  and  can  be  skipped.  Then, 
for  r  =  1  sin  bt, 


[(GiM2)T/2RT/2]T 


(  (  1  -  e~T/2) _ (z  sin  bT/2) _ 

( (z  -  e“T/2]  [(z  -  cos  bT/2)2  +  (sin  bT/2)2] 


Expanding  in  partial  fractions,  Eq.  A- 20  becomes 


(A-20) 


((  1  -  e~T/2  )  _ (z  si  bT/2) _ ) 

j(z  -  e-T/2)  [(z  -  COS  bT/2)2  +  (sin  bT/2 )2 ( 


Kiz  s in  bT/2 

A - -  +  - 

bT/2) 2+( sin  bT/2)2]  [(z  - 


cos 


K2z(z  -  cos  bT/2)  K3Z  ) 

cos  bT/2)2+(sin  bT/2)2]  (z-e~T/2)  } 


Solving  for  Kj ,  ,  and  in  Eq.  A-21  gives  Eq .  A-22  (z  -  es^2). 


(  (  L  -  e~T/2)  _ (z  sin  bT/2) _ 

|(z  -  e-T/2)  [(z  -  cos  bT/2)2  +  (sin  bT/2)2] 

j_  ii  -  g-T/2) _ 

_  cos  bT/2)2  +  (sin  bT/2)2] 


x  ( cos  bT/2  -  e  T/2jZ  sin  bT/2 
[(z  -  cos  bT/2)2  +  (sin  bT/2)2] 

_  (sin  bT/2)z(z  -  cos  bT/2)  (sin  bT/2)z  )  l 

[(z  -  cos  bT/2)2  +  (sin  bT/2)2]  (z  -  e~T/2)  j 


Since  we  are  concerned  with  steady-state  frequency  response,  only  the 
sine  and  cosine  terms  need  to  be  considered.  Therefore,  resampling  at  a 
T  frame  time  results  in  (z  »  esT): 


(  (  1— e~^-/ 2 )  (z  sin  bT/2)  j”^  _  ^lz  s^-n  ^T  Kyzfs  -  cos  bl ) 

((z-e~T/2J  [  (z  -  cos  bT/2)2+(sin  bT/2)2]j  [ (z  -  cos  bT)2  +  (sin  bT)2] 


where 


Ki  = 


(  .  .  . 

.12rU3 


’  /?  •>  c . 


—  T*  /  T> 
^  *■  i  j 


(e  -  cos  bT/2}^  +  (sin  bT/2)^ 


(A-23  ) 

(A-24  ) 


and 


k2 


(1  -  e~T/2)(sln  bT/2) 

"(e-^/2  -  cos  bT/2)2  +  (sin  bT/2)2 


(A-25) 


It  is  now  apparent  that  the  [  (G  if^)^/2]  ^  term  represents  a  sinusoidal 
forcing  function  for  Eq.  A-6.  In  Ref.  1,  Appendix  E,  it  is  shown  that 
Eq.  A--23  reduces  to  a  complex  multiplication  factor  in  the  frequency 
response  expression.  That  is,  the  complete  frequency  response  expres¬ 
sion  for  Eq.  A-6  is  given  by 


A-6 


A„  -I  j  11, 


A„  -1  j  B„ 

c  1  m2 
1/2 

| 

<h"2 

T 

<c2m)t/2 

,SPsT/2 

ji  +  1(g1m2)(g2M}t721 

rl-1 

1 

7,*CS^' 

[Kj + jk2] 

1  -■== 

7-1  oik(T/2) 

T 

(A-26) 

whe  re 

UJn  = 

.  2ti  n 

b  +  f72 

> 

n  0  j  1)  2 )  *  •  • 

(A-27 ) 

“k  = 

b  ■+ 

2nk 

T 

> 

k  -  0,  1 ,  2,  ... 

The  complex 

mult i p licat Ion 

factor  [K|  +  jK2 I  is  defined 

in. 

Eqs . 

A- 2  4  and 

A-  2  5 . 

The 

remai n 

fug 

f  requency 

response  factors  were 

written 

directly 

from  inspection  of  the  output  equation  (Eq.  A-6)  using  the  fre¬ 

quency  response  identities  presented  in  Section  IV,  Subsection  D. 

Notice  that  Kj  and  K.-,  are  functions  of  the  input  frequency  b  and  not 
=  b  +  2irk/T.  In  addition,  both  Kj  and  K2  are  periodic  functions  at 
b  =  ;'r/(T/2),  the  input  sampling  frequency.  As  discussed  in  Section  IV, 
the  result  of  this  singular  dependence  on  the  input  excitation  frequency 
b  and  the  periodic  nature  of  the  and  K2  functions  is  multiple  Bode 
plots. 

Substituting  the  previously  calculated  terms  into  Eq.  A-26  results 
in  Eq.  A-28: 


An+jBn 


( i-e-sr/2 ) 

T/2s(s+l) 


(l-e-sT/2) 

(  1-e  1 / 2 )  z+1 

Ts (s+1 ) 

s  =  J<»>n 

s  =  jwk  (z-e~T/2)  z 

“  '  - 

z=14u)k(T/2) 


(z  -  0.81)2 _ 

z2  -  1 . 61 z  +  0.6822 


z#esT 

Ul+jK2] 

z=  l4o)kT 


(A-28) 


This  equation  gives  exactly  the  same  results  as  Eqs.  159  or  175  in  Sec¬ 


A-7 


tion  IV. 


FREQUENCY  RESPONSE  TERM  [GRT/2]T 

It  Is  possible  to  derive  a  general  frequency  response  identity 
factor  for  the  term  in  Eq.  A-6.  Consider  the  more  general 
expression  in  Eq.  A-29  with  r  =  1  sin  bt. 


[GrT/M]T  =  {cT/MRT/MjT  (A-29) 


The  T/M  transform  for  the  sinusoidal  input  is  expressed  as 


rT/m 


•z  sin  bT/M 
-  (cos  bT/M)z  +  1 


z  =  esT/M  (A-30) 


Substituting  Eq.  A-30  into  Eq.  A-29  and  expanding  the  T/M  transform 
product  Into  partial  fractions  gives 


gt/m 


z  sin  bT/M 

T 

K^Z  sin  bT/M  +  K2z(z  -  cos  bT/M)' 

z^  -  (cos  bT/M)z  +  1 

z2  -  (cos  bT/M)z  +  1 

+ 


Terras  due  to 


T 


modes  of  G^/M 


(A- 31 ) 


Since  we  are  concerned  with  the  steady-state  frequency  response,  the 
partial  fraction  expansion  terms  due  to  the  modes  of  G^1,1  can  be 
dropped.  Therefore,  multiply  both  sides  of  Eq.  A-31  by  z2+(cos  bT/M)z+l 
and  evaluate  at  z  =  cos  bT/M  +  j  sin  bT/M.  The  result  is 


appendix  b 


COMPUTATION  OF  (W^M^)1  AND  (w2*M3W3)T 

As  noted  in  Section  IV,  the  computation  of  terms  such  as  (W3AM2W2)T 
and  (W2*M3W3)t  is  straightforward.  First,  evaluate  (W3*M2W2)T  using 
Eq.  210: 


(B-3) 


B-l 


Clearing  through  and  writing  Eq.  B-3  in  terms  of  the  advanced  2-transform 
for  1/s: 


-sT  s2T/3 

|  1  e  e 

-sT  s2T/3 
e  e 

-sT  sT/3 
e  e 

-sT  sT/3  ~sT 

e  e  e 

s 

s 

s 

s 

s  s 

sT/2 

e 

sT/6 

e 

sT/6 

e 

~sT  s5T/6 
e  e 

sT  s5T/6  -sT  sT/ 2 

e  e  e  e 

s 

s 

s 

s 

s  s 

1  0  0 

0  1  0 


(B-34) 


Equations  B-2  and  B-4  justify  the  assignments  made  in  Eq.  211  of  Sec¬ 
tion  IV. 


B-2 


I 


r 


APPENDIX  C 

NUMERICAL  DETAILS  FOR  AN  ILLUSTRATIVE  EXAMPLE 

Equation  223  of  Section  IV  stated 


vr/3 


- 


p^z 


,  ~T/3)  ,  -T/6.  ,  2  ^  -T/6  ,  -2T/6.  .  2, 

(p  -  e  <p  -  e  )(p  +  e  p  +  e  Hz  -  p  ) 


2,2  -5T/6) 

z  (z  +  e _ 


,  -2T/3 .  ,  -T/3.  ■  2  ,  -T/3  J  -2T/3. 

(z-e  )  (z  -  e  )u  +  e  z  +  e  ) 


z  = 

(C-l) 


There  are  four  residues  required: 


3 

_P_5_ 


|p-e 


|p-e 


,  -T  -T/3.,  -2T/3. 

Ce  -  e  )(z  -  e  ) 


-T/2 
e  2. 


(e-T/S  -  e-T/3)(3e'T/3)(z  -  e'T'3) 


(C-2) 


(C— 3) 


res,  --  res 


p—  (o  T/6/2)  +  (j/3/2)e 


-T/6 


K-  *  Lf  •-T,b)’ 


(-c~Tfb  ,  j/J  -T/6 


-  e-T/3)(^  +  If  cT1'*  -  e-T/6)(3/3  e'^6)  [z  -  (^-  +  ^  e-T'6) 


(C-4) 


C-l 


Switch  decomposition  can  be  used  as  an  alternative  (refer  to  Fig.  C-l) 


(T 

W3*  7+2  W2 )  (W2i 


R) 


=  w., 


( 

1 

y 

rt 

] 

sT/3 

1 

1  i 

\ 

e 

s  +  2 

( 

s2T/3 

e 

( 

_(cST/2R)T_ 

(C-5) 


=  W, 


Z 

-T 

e 

2  .  -3T/2 

2  4-  e  z 

-2T/3 

ze 

-5T/3 

e 

z 

-2T/3  2  ,  -13T/16 

e  z  +  e 

-4  T/3 

L  ze 

-T/3 
ze  J 

-T 

z.  -  e 

=  W3 

2,  -4 T/3  ,  -5T/6 

L  z  (e  +  e  J 

z  -  e 

-2T 

-T/2 

ze 

-2T.  .  -T. 

U  -  e  Hz  -  e  ) 

-T 

z  -  e 

(C 

Figure  C-l.  Switch  Decomposition  of  T/2,  T/3  Example 

Next,  define  z  =  and  clear  through  by  W_ : 


f. 


6  ,  -3T/2  3 

z  +  e  z 


6,  -2T/3.  ,  -13T/16  3 


m  /  O  /  T  /  Q 

Expand  Eq,  C-6  and  cancel  the  pole/zero  pair  z'-  +  e  '  +  e  : 


6  -2T/3  5  .  -4T/3  -ST/6  4  -3T/2  3  ,  -13T/1 6  2 

z  +  e _ z  +  (e  +  e  )  z  +  e  z+e  z 

,  -2.T/2,  ,  2  -2T/3  ,  -4T/3.  ,  -T/3w  2  ^  -T/3  -2T/3, 

(z-e  )(z+e  z  +  e  )(z-e  )(z+e  z+e  ) 


2,2  -5T/ 6. 

z  (z  +  e _ ) 


-2T/3.  ,  -T/3.  ,  2  -T/3  -2T/3. 

-e  )(z-e  )(z+e  z  +  e  ) 


C-3 


APPENDIX  D 


LISTING  OF  TRANSFER  FUNCTIONS 


The  s-  and  w'  transfer  functions  are  listed  for  a  T  =  0.04  seconds 
(25  Hz).  The  actuator  model,  20/(s  +  20),  is  included.  The  w'  transfer 
functions  presume  a  ZOH  and  were  calculated  using  the  DISCRET  computer 
program  presented  in  Volume  II. 

Also  listed  are  the  advanced  w'-transfer  functions  for  a  needed  for 
the  analyses  of  Section  V.  The  denominator,  which  is  the  same  regardless 
of  the  amount  of  time  advance,  is  listed  only  once.  The  numerators  are 
listed  for  advances  from  0.0  to  0.9. 

For  ease  of  presentation,  the  variable  s  has  been  used  in  all  of  the 
tables.  The  particular  data  heading  makes  clear  whether  s  or  w'  is  in¬ 
tended. 

s-  AND  w '-PLANE  TRANSFER  FUNCTIONS 

a/6e,  s- Plane 

AK=  (  ank/adk)  =  -1  . 44-9 
ATF(S)  NUMERATOR 


I 

AN POLY ( T ) 

AZKROf  T) 

1 

( 

-1 . 449 

)s**  5 

( 

-.8005E-02) 

+  J(  -.7455E-01) 

2 

( 

-1  47.6 

)  S**  2 

( 

-.8005E-0?) 

+  J{  • 7455R-01 ) 

2 

( 

—2.570 

)s**  1 

/ 

-101.9  ) 

+  J(  0.  ) 

4 

( 

-.8295 

) 

ANK= 

-1 .449 

ATF(S)  DENOMINATOR 


I  ADPOLY(T)  APOLE(I) 


1 

( 

1  .000 

)  S**  9 

(  1.157  ) 

+  J( 

0.  ) 

2 

( 

20.95 

)S**  4 

(  -.9545K-02) 

+  J( 

, 7705E-01 ) 

5 

( 

10.71 

)S**  5 

(  -.9545E-0?) 

+  J( 

— .7705E— 01 ) 

4 

( 

-bC.GC 

)S**  2 

(  -2.071  ) 

+  J( 

0.  ) 

9 

( 

-.7820 

)S**  1 

(  -20.00  ) 

+  J( 

0.  ) 

U 

( 

-.2857 

) 

ADK= 

1  . 

,000 

D-l 


a/8e,  »' -Plane 


RK=  (  RNK'RDK  )=  -.2926E-03 

RTF(S)  NUMERATOR 

I  RNPOLY ( I )  RZERO(I) 

1  (  — .  2926E-03  )  S**  5  (  -.8003E-02)  +  J(  .7455E-01) 

2  (  .I+986B-O1 )  S**  4  (  -.8003E-02)  +  J(  — .7455E-OI) 

5  (  1.040  )S**  3  (  50.00  )  +  j(  .1707E-10) 

4  (  — l40.0  )S**  2  (  -54.67  )  +  j(  0.  ) 

5  (  -2.236  )S**  1  (  175-1  )  +  J(  0.  ) 

6  (  -.7873  ) 


RTF(S)  DENOMINATOR 

I  KDPOLY(l)  -  -  -  ..  RPOLE(I) 

1  (  1 .000  )S**  5  (  1-137  )  +  J(  0.  ) 

2  (  19-95  )S**  h  (  -.9345E-O2)  +  J(  .7705&-O1) 

3  (  15.74  )S**  3  (  -.9345E-02)  +  j(  -.7705E-01) 

4  (  -44.29  )S**  7  (  2.070  )  +  ,J(  0.  ) 

5  (  -.7430  ) S**  1  (  -19.00  )  +  J(  0.  ) 

6  (  -.2693  ) 


q/6e,  s- Plane 

BK-  (  BNK'BDK  )=  — 1 47.6 

BTF(S)  NUMERATOR 

I  BNPOLY ( I )  BZERO(I) 

1  (  -147.6  )S**  3(0.  )  +  J(  0.  ) 

2  (  -79-39  2  (  -.88160-0?)  +  J(  0.  ) 

3  (  -.6884  )S**  1  (  -.5292  _  )  +  J(  0.  ) 

4  (  v.'.  )  BNK-  — 147.6 


BTF(S)  DENOMINATOR 

I  BDPOLY(l)  BPOLE(l) 

1  (  1 .000  )S**  5  (  1.137  )  +  J(  0.  ) 

2  (  20.95  )S**  4  (  -.9345E-02)  +  J (  .7705E-O1) 

5  (  16.71  ) S**  3  (  -.9345E-02)  +  J(  -.7705E-01) 

if  (  -46.66  )s**  2  (  -2.071  )  +  J(  0.  ) 

5  (  -.-7820  )s**  1  (  -20.00  )  +  j(  0.  ) 

6  (  -.2837  )  DDK-  1  .COO 


R-2 


q/5ft,  v' -Plane 

3K=  (  SNK'SDK  )=  .Y546E-02 

oTF( S)  NUMERATOR 

]  snpoly(i)  szero(i) 

1  (  .7546K-02)S**  5  (  —.88l  6E— 02)  +  j(  0.  ) 

2  (  ?.4p8  )S*»  4  (  --529?  )  +  J(  0. 

3  (  -138.8  ) S**  3  (  50.00  )  +  J(  -.1717E-11) 

J,  (  -73.34  )s**  ?  (  -571.2  )  +  J(  0.  ) 

5  (  -,6534  )  S**  l(0.  )  +  J(  0.  ) 

C  (  0.  ) 


3TF(  S')  DENOMINATOR 

1  SDfOLY(l')  SPOLK(I) 

1  (  1.000  ) 3**  5  (  1-137  >  +  l 

p  (  19.Q1'  is**  4  (  —  .9343E-02)  +  J(  . 7705B-01 ) 

3  (  T3.74  )s**  3  (  -.9345E-02)  +  J(  -.7705E-01 ) 

.4  (  -44.29  ) S**  2  (  -2.070  )  +  j(  0.  ) 

3  (  -.7430  'IS**  1  (  -19.00  )  +  J(  0.  ) 

C  (  -.2693  ) 


An/8e,  8- Plww» 

CK=  (  CNK/CDK  )=  -.5315 

ctf(s)  numerator 

1  cnpoly(i)  czero(i) 

1  (  -.3315  )s**  4(0.  )  +  j(  0. 

P  (  -.331 4  )s**  3  (  .1831E-02)  +  J(  0. 

3  (  -24.95  is**  2  (  -.5009  )  +  j(  35 

4  (  .4567E-01 )s**  1  (  -.5009  )  +  0.853 

iz  (  n  CNK-  .5315 


) 

) 

> 

; 

) 


cte(s)  denominator 

1  CDTOLY( 1)  CPOLE(l) 

1  (  1.000  )S**  5  (  1.157  )  +  J(  0. 

(  po.95  )s**  4  (  -.934551-02)  +  J(  .7705E-01 

3  (  16.71  )s**  3  (  -.9345E-0?)  +  J(  — 7705E-01) 

4  (  -46.66  )s**  2  (  -2.071  )  +  J(  0.  ) 

3  (  -.7820  is**  1  (  -20.00  )  +  J(  O.  ) 

6  (  —.2837  i  CDK=  1.000 


D-3 


t 


An/&e;  V -Plane 


TK=  (  TNK/'TDK  )=  .  1 003E-O1 

ttf(s)  numerator 


I 

tnpoly(i) 

TZERO( I) 

1 

( 

•1003E-0i)S**  5 

( 

.1830E-02)  +  j(  0. 

1 

2 

( 

-.490?  )S**  4 

( 

-.5721 

)  +  j(  6.847 

) 

3 

( 

-•9944e-Oi ) S**  3 

( 

--.5721 

)  +  J(  -6.847 

) 

4 

( 

-23.68  )S**  2 

( 

50,00 

)  +  J(  — . 2797E-1 1 ) 

5 

( 

•4333b-oi)s**  1 

( 

0. 

)  +  J(  0. 

) 

6 

( 

0.  ) 

! 


ttf(s)  denominator 


tdpoly( I) 


TPOLE( I) 


1 

( 

i  .000 

)S**  5  ( 

1.137  ') 

+ 

J( 

0. 

)  ■  . 

2 

( 

19.95 

) S**  4  ( 

-.9345E-02) 

J( 

.7705E-01  ) 

3 

( 

15.74 

) S**  3  ( 

-.9345B-02) 

+ 

<T( 

-.7705E-01) 

4 

( 

-44,29 

)  S~*  2  ( 

-2.070  ) 

+ 

J( 

0. 

5 

( 

-.7430 

)S**  1  ( 

-19.00  ) 

4* 

J( 

c. 

' 

6 

( 

-,2693 

) 

1 

a  TIME  ADVANCED  TRANSFER  FUNCTIONS  ( w’ ) 

DENOMINATOR 


1 
1 

2 
3 
it 

5 

6 


npoly(i) 

(  —  .2926Fr-03)S**  5 

(  ,U986b-oi)s**  U 


(  - 
(  - 

(  - 


i  .o4o 
i4o.o 
2.  236 
7873 


)s**  3 

)S**  2 
)S**  1 
) 


NUMERATOR 

-.BOOJE-O?) 
-.8003E-02) 
50.00  ) 

—5^.67  ) 

175.1  ) 


D-4 


ZERO ( I) 


J( 

j( 

J( 

J( 

J( 


.7455^01) 
-.7455^-01 ) 
.1707E-I0) 
0.  ) 

0.  ) 


I 

DPOLY(I) 

POLE( I) 

1 

( 

1 .000 

)s**  5 

(  1.137  )  +  j(  0.  ) 

2 

( 

19.95 

)S**  4 

(  -.9345E-02)  +  J(  .7705E-01) 

*  | 

3 

( 

15.74 

)  s**'  5 

(  “.9345E-0?)  +  J(  --.7703E-01 ) 

4 

( 

-44.29 

)S**  2 

(  -2.070  )  +  J(  0.  ) 

5 

( 

-.7430 

)s**  1 

(  -19.00  )  +  j(  0,  ) 

6 

( 

-.2693 

) 

A  = 

0  (0.00  seconds) 

1 

A  =  .1  (.004  seconds) 
NUMERATOR 


1 
1 

2 
3 

it 

5 

6 


npoly(i) 

(  -.i474u-o3)s**  5 
(  .5291E-01  )S**  4 

(  .4797  )S**  5 

(  -i4o.o  )S**  2 

(  -2.239  )s**  i 

(  -.7873  ) 


ZERO ( I ) 

(  -.8003&-02)  +  J(  .7455E-01) 
(  -.80030-02)  +  j(  -.7455&-01) 

(  50.48  )  +  J(  0.  ) 

(  -32.18  )  +  J(  0.  ) 

(  360.7  )  +  j(  0.  ) 


A  =  .2  (.008  seconds) 
NUMERATOR 


I  NFOLY(l)  ZERO(l) 

1  (  .85250-07')  s**  5  (  -.80030-02)  +  J(  -.7455E-oi) 

2  {  .53720-01  )S>#  4  (  -.80030-02)  +  J(  .7455E-01) 

3  (  60560-01 )s**  3  (  51*02  )  +  J(  0.  ) 

4  (  -1 40.1  )S**  2  (  -50.31  )  +  j(  0.  ) 

5  (  -2.24-2  )S**  1  (  -.6301E+06)  +  j(  0.  ) 

6  (  -.7873  ) 


A  =  .3  (.012  seconds) 
NUMERATOR 


1  NPOLY(l)  ZERO(l) 

1  (  .1  4o6e-03)s**  5  (  -.60030-02)  +  J(  .74550-01 ) 

2  (  .52280-01 ) S**  4  (  -.80030-02)  +  J(  -.74350-01) 

3  (  -.64 08  )s**  3  (  -48.93  )  +  J(  0.  ) 

)i  /  _ 1* r\  «  S  n  (  c;)i  AO  ^  4.  Tf  Q  ^ 

\  I  ‘T'w  •  |  ;  U  •’  -  C.  \  9  wc-  /  7  w  \  V  -  / 

5  (  -2.245  )S**  1  (  -376.9  )  +  j(  0.  ) 

(  -.7873 


1 
1 

2 

3 

4 

5 

6 


A  =  .4  (.01 6  seconds) 


P0LY( I) 

(  .2656EM)3)S**  5 
(  . 48590-01 } S**  4 
(  -1.201  )S**  3 

(  — 1 40.1  )S**  2 

(  -2.249  )S**  1 

(  -.7873  ) 


ZER0(I) 

(  -.80030-02)  +  J( 

(  -.80030-02)  +  J( 

(  -47.95  )  +  J( 

(  57.22  )  +  J( 

(  -192.2  )  +  J( 


. 7455E-OI ) 
-.74550-01) 
0.  ) 
0.  ) 
0.  ) 


D-5 


cwji  -fv>i  rv>  -1  m  ct\\ji  -f-  ^  ro  -1  h  °S'VJ1 


A  =  ,‘j  (.02  seconds) 


NUMERATOR 

I  IMPOLY(l)  ZERO(l) 

1  (  .36722-03) S**  5  (  -.80032-02)  +  J(  .7455E-0 1) 

2  (  .  42662-01  )S**  4  (  -.800513-02)  +  J(  -.7455B-01) 

3  (  -1.761  )S**  3  (  -47.34  )  +  J(  0.  ) 

4  (  -140.1  )S**  2  (  61.72  )  +  J(  0.  ) 

(  -2.252  )S**  1  (  -130.5  )  +  J(  0.  ) 

(  -.7873  ) 


A  =  .6  (.024  seconds) 

NUMERATOR 

NPOLY(T)  zkro(i) 

(  .l+303E-O3)S**  5  (  -.80032-02)  +  J(  .74-55B-01) 
(  . 34482-01 ) S**  4  (  -.80032-02)  +  J(  -.74552-01) 
(  -2.322  )S**  3  (  -47.10  )  +  J(  0.  ) 

(  — 1 40.1  )S**  2  (  68.09  )  +  J(  0.  ) 

(  -2.235  js**  1  (  -99.64  )  +  J(  0.  ) 

(  -.7875  ) 


A  =  .7  (.028  seconds) 

NUMERATOR 

NFOLY( I)  ZERO( I) 

(  .47222-03) s**  5  (  -.80032-02)  +  j(  .74550-01 ) 

/  I .  1  1.  /  Qr\s\'2T'  r\r->\  .  T  (  r~7)i  — r>1  ^ 

V  •  ( rr-u  1  ;  ^  4  v  -.uw^u-vc/  -f  u\  •  1  ^  »  / 

(  -2.882  )S**  3  (  -47.23  )  +  j(  0.  ) 

(  — 1 4o. 1  )S**  2  (  77.41  )  +  J(  0.  ) 

(  -2.258  )S**  1  (  -8l.l3  )  +  J(  0.  ) 

(  -.7873  ) 


A  =  .8  (.032  seconds) 

NUMERATOR 

1  NPOLY(I)  ZERO(I) 

1  (  ,46?92-03)S**  5  (  -.80032-02)  +  J(  .74552-01) 

2  (  .1 1  4l  13-01  )S**  4  (  -.80032-02)  +  J(  -.74552-01) 

3  (  -3.442  )s**  3  (  -47.79  )  +  j(  0.  ) 

4  (  — 1 40- 1  )S**  2  (  -68.83  )  +  J(  0.  ) 

5  (  -2.261  )S**  1  (  91.98  )  +  j(  0.  ) 

G  (  -.7873  ) 
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A  -  .9  (.056  seconds) 


numerator 


1 
1 

2 

3 

4 


NPOLY { T) 

(  .4o47E-03)s**  5 

(  -.348lE-02)s**  4 
(  -4.003  )s**  3 

(  -l4o.i  )s**  2 

(  —2. ?64  ) 0**  1 

(  -.7873  ) 


\ 

( 

( 

( 

( 


ZERO(l) 

— .8o03B—02)  +  j( 

— .8oo3b-0?)  +  j( 
-j*8.87  )  +  j( 

-60.19  )  +  j( 

117-7  )  +  J( 


-7435E-01) 
-.7435B-01 ) 
0.  ) 
0.  ) 
0.  ) 


i 


D-7 
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